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ABSTRACT 


The spectral analysis technique is used in order to predict the 
wave induced rolling motions of a ship in a seaway. Non-linear 
equations of motion for the ship, the anti-rolling tank, and 
for wave excitations are developed. Solution techniques for 
the equations are examined, and a linearized form is adopted. 

A design method using a digital computer program based upon 
numerical solution of the linearized equations is proposed and 
examined for a number of different representative ships. The 
design method is demonstrated to be capable for use to design 
passive anti-rolling tanks for ships and to determine feasi- 


bility of activated anti-rolling tanks. 
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Nomenclature 


Area of fluid surface in one wing tank 
Sectional area of cross duct 


Area of any tank section perpendicular to l, 
a= A at 1-70, L 


Location of the center of buoyancy 
Beam of ship 


Distance from the center of buoyancy to the 
center of gravity 


Tank linear damping coefficient 
Wing tank breadth 
Boy (L B T) 


Block coefficient, 


Longitudinal Prismatic coefficient, Cb/ 
Midships sectional area coefficient 


Moment arm of section tank area a about G 
ea )/dt 

al )dt 

Sway damping ordinates 
Surface elevation of the sea, positive downwards 
Tank movement 

Tank excitations 

Sway damping frequency abscissas 


Location of the center of gravity of the ship 


Controlled feedback gains 


. Metacentric height (transverse) of the ship, 


corrected for all free surface except for that of 


the tank: 
Height of the wing tank 


Height of the cross duct 
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KB 


KG 


Le 


Lo 


Significant wave height, the average of the 1/3 
highest waves 


General subscript 

Squareroot oft -1 

Moment of inertia of the ship about the x-axis 
Wave constant, 63" /ሮ 

Roll moment 


First order roll velocity coefficient 


First order static righting coefficient 
First order roll acceleration coefficient 
Second order roll velocity coefficient 
Dimensionless ship damping in roll 
Dimensionless tank damping 

Dimensionless ship quadratic damping in roll 


Dimensionless cubic spring constant for ship in 
roll 


Dimensionless ship damping in sway 
Distance from the keel to the center of buoyancy 
Distance from the keel to the center of gravity 


Centerline distance along tank sections in the 
+n direction 


Endpoint of 1 
Length of ship 
Effective length of ship 


The maximum fluid excursion along 1 from 1 = 0 
in the Tr direction 


Wing tank length 


Origin of nmn coordinate 


vii 





p(t) 


p* 


Mass of ship 

Angular roll velocity, ¢ 

Angular roll acceleration, <? 

Pump pressure 

Dimensionless pump -pressure, p(t)/ p gR 
Radial distance from G to a point along 1 


Distance from the centerline of the ship to the 
center of area As 


Distance along wave from reference 
Weighted tank section area, dé [a 
O a 

Weighted tank moment arm, > [gjar 
O R 

Wave elevation spectral density function 

Time 

Pump lag time 

Draft to the center of buoyancy 

Draft of ship 

Speed of ship in the x-direction 

Sway velocity of ship 

Dimensionless sway velocity, v/wsB 

Horizontal wave velocity 

Tank weight 


Ship reference coordinate axes 


' Barth referenced coordinate axes 


Heading angle parameter, Le |cos (8) | /A 
Sway direction 


Sway force 


viii 





ei 


Bi 


Oi 


Yi 
Ab 
Ào 
Ày 
Nt 


First order sway damping coefficient 
First order sway inertia coefficient 
Distance from G to B in the positive z direction 


Sway coefficients in the expanded equations of 
motion 


Distance from G to 0 in positive z direction 


Parameter of the Pierson-Moskowitz Spectral Densi- 
ty Une Enron 


Dimensionless coupling parameters in the expanded 
equations of motion 


Parameter of the Pierson-Moskowitz Spectral 
Density function 


Heading angle, 180° for head seas 


Dimensionless coefficients of roll in the expanded 
equations of motion 


Dimensionless tank coefficients in the expanded 
equations of motion 


Displacement of ship, long tons (2240 lb) 
Increment of the parameter 


1/3 ( tank)/ 1 


Tank effectiveness parameter, 6 
max max tank) 


Roll angle from the vertical 
Maximum static equilibrium angle 
Waveslope angle 

Inertia constant 

Ship-sway coupling parameter 

Va Sway coupling parameter 

V. Sway coupling parameter 


Ship-tank coupling parameter 
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Dimensionless driving frequency, W/W 
Tank tuning parameter, OW, / 0. 


Ship-tank decoupling frequency, 034/02 
Sway decoupling frequency, Y9/B/0. 
Frequency 


Encounter frequency of the waves as seen by 
the moving ship 


Absolute frequency relative to the earth 
sShıpsnatural GUN frequency, VK9/(K, - Ix) 
Ship-tank decoupling frequency 

Tank natural frequency 

Yg/B 

Mass density of tank fluid 

Mass density of seawater, 35 ር” 


Dimensionless time, wt 


Note: Computer program symbols are given 
separately in Appendix C 





fen introduction 
1. General Discussion 
Successful roll stabilization of ships by means of 
fluid filled Anti-Rolling tanks requires a thorough analysis 
of the forces and moments involved. The author undertook 
this task with a goal of being able to use the results as a 
design tool of the naval architect. As such, the Anti-Rolling 
tank design method derived in this paper is intended to be 
a single step in the larger process of designing the whole 
ship. 
Modeling the dynamics of a ship with an Anti-Rolling 
tank has been investigated by many authors of whom only a 
small percentage are represented in the reference list.  Sev- 
eral problems have prevented their analyses from being complete 
or from possessing sufficient accuracy in the range of interest. 
One central problem is that rolling is a resonance 
phenomenon. The ship is an underdamped dynamic system in roll, 
one which magnifies wave excitations to a certain degree. This 
may be proved by considering that it is well established for 
waves in the open ocean to have a maximum waveslope of about 
l6 degrees. Waves much steeper than 16 degrees are unstable 
and break up by forming whitecaps or similar turbulence. 
However, rolls of 20 or 30 degress or more are common even in 
moderate seas. Thus, for the normal range of ship parameters 
and dimensions, a resonant condition occurs within the fre- 
quency range of wave exciting forces. 


To better understand the implications of being ina 





resonant condition, consider the ship as a simple second 
order system in roll, with the rin ጆች 0 representing 
roll angie. 

(Kp-Ix)"$ + Kpó + Kod = - (Wave excitations) 


The natural resonant frequency of the system is w_= Ko 
i | Upon 


and it may be recognized that the amplitude of the roll motion 
at the resonant frequency is proportional to VKd (Kp-Ix)/Kp 
Hence the magnitude of the rolling motion at resonance will be 
most strongly dependent upon the damping parameter Kp. Some 
analyses of ship motions have assumed a linear damping para- 
meter for the ship, i.e. Kp = a constant. This is a good 
approximation for small motions or for specific conditions, 
but in the general case we expect Kp to be non-linear. Since 
roll damping involves viscous friction, we expect that there 
will be a Reynold’ number dependency. Also, because rolling 
creates waves, there is added damping due to the radiation of 
E disturbances from the ship. This radiation of wave energy 
has been shown to be frequency dependent in a complicated 
fashion. 

Enocher problem is that rolling motion is not indepen- 
dent of other ship motions and therefore cannot be reasonably 
expressed as a simple second order system. For example, a 
ship at rest has been caused to assume an angle of roll which 
is different from its usual upright position. This roll angle 
causes the ship to have an unsymmetrical underwater shape, 
which induces a righting moment (assuming a stable ship). 


Obviously, any pitching motion would change the underwater shape 





still further, unless the ship were prefectly spherical, 

and would induce a change in the value of the righting moment. 
It can be shown by similar arguments that roll is coupled to 
all of the other ship motions and displacements. However, 
extensive studies have shown that the degree of coupling is 
negligible in all instances, except for sway rate and sway 
acceleration, assuming a "normal" surface ship. 

This paper will take an engineering approach towards 
analyzing rolling motion. The best available information on 
how to model non-linearities will be used in a workable 
system which allows calculations to be done in a reasonable 
amount of time. This approach was taken because of the 
iterative nature of the design process for Anti-Rolling tanks, 
which in turn is one step in a larger iterative design process 
for the entire ship. Therefore the analytical model of the 
Anti-Rolling tank had to be flexible enough to be "turned 
inside out" into a systhesis model by which a large range of 
parameters could be investigated by some search technique. 
fais factor was the primary motivation behind taking an engin- 
eering approach to a sufficient design rather than a purely 
analytical approach towards a perfect dynamic model. It is 
hoped that this paper possesses a good enough balance between 
numerical estimation and reality to shed a little more light 
on a very complicated subject. 

2. System of Units 
The English system of units is used throughout this 


paper and in the computer program. The following English units 





are assumed pertinent to all applicable quantities unless 
other units are indicated. 

Time: seconds 

Distance: feet 

Mass: slugs 

Angle: radians 

Gravity constant: 32.2 ft. per sec. 

Displacement: long tons (2240 lbs) 

Force: pounds 


Moment:  foot-pounds 





TE. Derivation of the Equations of Motion 
1. General Discussion 

This chapter presents a detailed derivation of the 
Pemations of motion of the ship-tank system. The form of 
the equations is highly dependent upon the coordinate system 
chosen, which makes it important to understand section 2 and 
figures l and 2 before continuing with the text. For those 
who are interested only in the results, section 9 presents 
the final form of the equations and section 10 gives their 
non-dimensionalized equivalent. It is suggested that readers 
who are not interested in the derivations should still read 
section 2 before skipping to the results. This will ensure 
the reader has a reasonable feeling for the basic nomenclature 
of the equations before facing them. 

The derivation of the equations of motion begins with 
Newton's second law of motion, action equals reaction. All 
of the terms which make up the reactive forces and moments are 
first listed. These terms are each broken down into their 
dependent variables. An expression for the reactions is then 
formed by taking the first partial derivative of the partic- 
ular force or moment with respect to each of it's dependent 
variables. If the indicated partial derivatives were all 
ሺ ከ. nts the reaction terms would be exactly linear. Each of 
the partial derivatives is examined for this property, simpli- 
fications are made where possible and in several cases secondary 
parameters are introduced to further define particular coeffi- 


cients. Enough documentation of relevant assumptions is made 





to be able tto change particular coefficient values whenever 
the state of the art permits. The exei ton forces and 
moments are handled in the same fashion except that they 
have independent as well as dependent variables. 

The last section of this chapter gives the summarized 
results as a set of simultaneous equations which are in non- 
dimensional form. This section is the key reference in 
understanding the rest of this paper. This is because the 
non-dimensional coefficient parameters found there constitute 
a kind of universal set which not only appears in most other 
literature on Anti-Rolling tanks (often with the same symbology) 
but also is used continuously throughout the rest of the paper. 
2. Coordinate System 

BtEacoobsdinate axes forzthe ship are through the center 
of gravity of the ship and lie along the principal axes of 
inertia. They are x,y, and Zz. The reference coordinates 
with respect to the earth are ee and z i positive z 


O 


erection is in the direction of gravity. 





2 


Figure 1. Ship-Oriented Coordinates 


The tank coordinate, n, is relative to the ship and is 


positive in the direction shown on Figure 2. All of the symbols 





shown on Figure 2 are in their positive sense. 

Duesgnteunteor this figure is to show the relative 
positions of points G, B, and O. G denotes the center of 
gravity of the ship with a filled Anti-Rolling tank (cross 
duct blocked) and serves as the reference for the ship 
oriented x-y-z axes. B, the center of buoyancy according to 
Archimede's principle, is assumed to be the point of appli- 
cation of the wave excitation forces and moments. O is the 
origin of the tank coordinate n. Also shown is the waveslope 
angle Y and some of the basic Anti-Rolling tank geometry. 
Figure A-1 in appendix A gives a more complete picture of the 


Anti-Rolling tank by itself. 
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3. Reaction Forces and Moments. 

For the ship-oriented coordinates a force and moment 
ር ርር 156 done. -K denotes a roll moment about the x-axis, 
Y denotes the force in the y direction, and E denotes the 


tank moment about an axis parallel to the x-axis. 


E (Moments in $ direction) = K = I? 


£ (Forces in the y direction) = Y = mv 


For the tank coordinate, reference (4) uses a Lagrangian 
description of the potential, kinetic, and dissipation 
energies associated with the tank to arrive at its equation 
of motion. For the sake of brevity, only the results, 
which are well known, will be used in this paper. 

Assuming a full tank with blocked ducts, there will 
be no interaction between the tank and roll or sway. In 
all other cases, there will be coupling between the three 


coordinates. This leads to an expression for K and Y as 


follows: 
BU -K(tank) + K(hydrodynamic) + Klexcitations) 
Y = Y (tank) + Y(hydrodynamic) + Y(excitations) 


The hydrodynamic forces and moments are those imparted 
to the body of the ship by virtue of the ship motions 
Occuring within and at the surface of a dense fluid (sea- 


water) in the presence of gravity. The tank forces and 
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moments are caused by the tank fluid shifting its weight 
transversely in the ship. Free surface waves at the surface 
of the ocean are the cause of excitation forces and moments. 
This separation of terms in the force and moment 
equations may be carried further by noting which variables 


correspond to tank, hydrodynamic or excitation portions. 


K or Y(tank) = K or Y (n, f, n) 
Kor Y(hydrodynamic) = K or Y (0, inp "n awe) 


ky V (exzciliations) = K or Y (9) 


The usual process of obtaining a first order model 
KE system is to take the first partial derivatives of each ` 
quantity (K or Y in this case) with respect to each of 
its arguments. After following the above procedure, 
combining terms where possible, and separating all the reactive 
forces and moments, including those of the tank, from the 


excitations, the roll and sway equations are obtained. 


(roll) (KyD + Kv)v + [(Kp - Ix)D? + KpD + Ko]@ + K(tank) = 
-K (excitations) 

(sway) [(Y% - m)D + Yv]v + YpD?2 + ypD + Y¢)¢ + Y(tank) = 
-Y (excitations) 


where: D = d/dt 


Kv = 3/(3v)K Yv = 9/(3v)Y 
Aou etc. 
Kp = 3/(8p)K 





Tl 


Kp 


9/ (95) 


Ko 


3/(39)K 


The first order equation of .motion for the Anti-Rolling 


tank, as found in reference (4) is: 


(tank) pA RDV - (S"D2 + Days - (sup =: B,D/2pA R“ ES 
29)01 55...) 
S' and S" are geometry factors which are calculated 


in Appendix A. 
The coupling terms between roll, sway, and the tank are: 


Knn + Knn + Knn 


K (tank) 


그 vom 


Since the coupling terms must be symmetric in any energy 
conservative system, the tank equation gives us directly 
the indicated partial derivatives in the roll and sway 


equations. 


Kn = -pA¿R?S" Yn = +2pA R? 
En = 0 Vins 
Kn = -2pgA_R? Yn = 0 


Using these results and taking the case of a purely 
passive Anti-Rolling tank (E. (t) = 0) the linear equations 


of motion for the ship-tank system are formed. The result is 
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a fifth order differential equation which is linear if the 
assumption is made that the COefficients of all the terms 


are constant. 


(roll) (KVD + Kv)v + [Kp - Ix) D? + KpD + Ko]¢ ~- oA, R2 (s"p* t 
ac 6 
° 2 ° 2 2 2 == 
Koa s [(Yv - m)/D + Yvjv + (YpD” + YpD)¢ + 2pA,R D°n = 
EXC) 


(tank) 2p AgGR?Dv ፦ pAGR2 (S"D* + 2000 =- pA R* (S'D2 + By / 


pAoRZD + 2g)n = -E,(t) = 0 


Ti order to add a control system to the tank (mathematically) 
it is necessary to set the right hand side of the tank 
equation, -Eg(t), equal to zAeRD(Ut); —p(t) 16 the pump 
pressure developed across a control effector in the tank 
duct, with a (+)p direction in the (-)n direction. A 
separate equation for the dynamics of the pump must be 
included in the equations of motion. A simplified form of 
p(t) which should be sufficient in the general case is to 
consider that any control effector can be characterized in 
a gross sense by an equivalent first-order time lag. This 


leads to the equation: 


T P(t) + p(t) = polt) 
Pe \t) 3 F (os es mi ni b) 


919 + g2% + gpn + 940 ፕ አ ን 그 gcV ug. m Pmax 


ㅣ 


pet) 


*Pmax: |919 *..-*9869| > Pmax 





ፐን 


new ቪን ን ሙቢ ያ. [Orpi (Mwasi hosen based on linear 
regulator theory which allows feedback of all the state 
variables, plus the excitation y. In practice, w is 
random and unmeasureable sona si Eny and of the state 
variables, the roll and tank variables are the easiest to 
measure and the most strongly coupled. This reduces the 


controller equation down to: 


p.(t) = gj$ + 996 + gan + 941፣ |8..| < Prax 


De) pr 1970 + god + gan + gan| > Pmax 


A useful equation which must be considered along with 
the above set of state equations is an estimation of the 
power requirements of the control effector. Using the 
relationship that power is the product of pressure differ- 
ential and volume flow, and taking into account pump efficiency, 


power requirements are calculated by: 


Pump Power (1/ey) p (t) AgRn 


where: eb - 0.5 - 0.8 


imp li fications of the Hydrodynamic Coefficients. 


Assume that the roll and sway motions of the ship 
as seen at the center of buoyance are uncoupled; that is, 
a roll motion about the center of buoyancy induces no 
Swaying and likewise a sway motion at the center of buoyancy 


induces no rolling there. It may be easily imagined that 
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this is nearly true. Using this basic assumption, many 
simplifications in the terms of the hydrodynamic coefficients 
are possible. For example, a roll motion about the x-axis, 
which passes through the center of gravity, will be seen at 
the center of buoyancy as a roll motion plus a sway motion. 
The hydrodynamic resistance to sway at the center of buoyancy 
is felt as a sway force, but does not couple back into roll, 
with this assumption. Usina the (T) notation to indicate 
that a particular coefficient is taken at the center of 
buoyancy instead of being referenced to the center of gravity 
(which is the center of the roll and sway axes) the values 


of the following hydrodynamic coefficients are determined. 


BEV — Yv 2 Oe 2 YY 
Y'v - vv Yp = Y'p = -z Yy 
K'd = Kb Kv = K'v = -ZQYv 

Kv = K'v = -2 Yv 


Introduction of the above relationships shows which 
hydrodynamic coefficients must be known as a minimum in 
order to solve the equations of motion of the ship-tank 
system. Notice that these coefficients are all defined 
as being determined with respect to the x-y-z coordinate 


System about the center of gravity of the ship. 


Kó Static righting moment characteristic of the ship. 


Kp Linear damping in roll. 
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Kp = Added inertia of the ship in roll to accomodate water 
flow. 


Yv = Linear damping in sway. 


Yv = Added inertia of the ship in sway to accomodate 
| water flow. 


5. lank Coefficients. 
The parameters which determine the coefficients of the 


tank variables in the linear equations of motion are: 


pA¿R? = Characteristic of tank size. 


OT = Function of tank geometry = 5? (a/R) al 
S = Function of tank geometry = £3 (Ag/a) al 
B, — Tank linear damping. 


S' and S" are the weighted moment arm and sectional area, 
respectively, of the tank fluid. They are familiar terms 
in literature on hydraulic engineering. oA, R? has the units 
of mass times a moment arm and is a characteristic tank 
size or inertia term. 

The above tank parameters are all free to be መን 
However, it can be noticed that S" and S' fall within the 


following range of possible values. 


= 
|^ 


S e 
-1 < S" < 1 


0.05 (V2957 pApR) 


|” 





16 


For the case of S' getting very large, the tank dynamics slow 
down due to the increased fluid inertia. The extreme case 
is the equivalent of the tank being just an additional 
Peronteaboard ship. The sign change of S" is accomplished 

by changing the center of gravity of the tank relative 
vertically to the center of gravity of the ship. A negative 
sign indicates a high tank, a zero value occurs near the 
center of gravity of the ship, and a positive sign is caused 
by a low tank. 

The damping parameter of the tank, 547 is free to be 
chosen as long as it exceeds the approximate minimum shown 
above, say about 5 percent of critical damping. Later 
analyses will show that a good passive tank will incorporate 
far more damping than this. Otherwise, undesirable secondary 
resonance peaks appear in the frequency response of the 
ship-tank system. Therefore, for practical purposes B, may 
be considered a free parameter with the stipulation that 
in the final design of the tank suitable damping devices such 
as sharp corners, baffles, etc. are used to achieve the 
desired degree of energy ጋ ር By may be calculated 
or predicted accurately enough with standard hydraulic 


engineering 'pipe friction! tables. 


6. Effector Pump Parameters. 


The maximum efficient effector size can be easily 





17 


estimated for any particular tank geometry, since the effector 
pump is installed in the cross ducting. “his will enable 
the choice of a particular sized pump which we have chosen 


to characterize by its maximum pressure, and lag time, 


Pmax’ 


E There is no particular requirement for p except 


max’ 
that 이 = (Power, max.)ep/AgRn max. For 22 however there 
is the requirement that the pump time lag must be less than 
me natural resonant frequency of the tank for the pump 


to have much useful effect. A pump with more lag will 


introduce instabilities into the control system. 


7. Non-Linearities. 

The hydrodynamic coefficients obviously all depend upon 
ship speed because they are caused by changing the flow 
field of the water about the ship. Of theSe, the damping 
-parameters, Yv and Kp, are the most speed dependent because 
they involve the dynamic lift in the horizontal plane generated 
by the ship hull serving as an inefficient lifting body. 
An angle of attack is imparted to the mean water flow by 
the relative motion between the ship and the water in the 
horizontal plane determined by the vector addition of sway 
rate and forward speed. The resultant lift force is felt 
as a resisting force in the sway direction. The static 
righting moment coefficient, K$, is the least speed dependent 
because the only change from zero speed conditions is that the 


constant water pressure surfaces, such as the air-water inter- 
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face, are no longer planar, but are slightly folded into a wave 
pattern. Based on these arguments, and with no further proof, 

it was decided to adopt Kp and Yv as ship speed dependent, and 

the other hydrodynamic coefficients as speed independent. 

Because the magnitude of the roll angle at resonant fre- 
quenty is much larger Se its overage yo lue i t helvalue of Kg 
does not remain reasonably constant, say within 5-10%, within 
its expected range. Since the well known curve of Static K 
versus $ is anti-symmetric, it was decided that it would be suf- 
ficient to adopt an odd polynomial expression for Kg, as was 
done in reference (7). The Kg% term is replaced by: 

Kg + 303 

L costficient aay be simply evaluated by noting that 
the value of K is zero when ø is at the maximum static equilibrium 
angle, ØR . This fixes the value of Kg3 as -Kø/ Øp? ThS can 
Be E uneturther szmplified by approximately gp by about 80° or 
radians. This enables the use of Kg3 — -Kg/2 if better in- 
formation is lacking. 

Part of the.energy dissipation caused by roll motion is 
Be) to viscous friction resisting the relative motion of the 
skin of the ship to the water. Since viscous friction is Rey- 
Beles Number dependent, it must be roll velocity dependent. This 
component of roll damping can be calculated by standard hydraulic 
መስክ ይፐ 'friction factor' tables to include skin friction 
and the eddy-making resistance of the bilge keels. However, 


this usage precludes us from making the odd polynomial 
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K T CIO because the standard tables are all of the form 
1/2 (area) (velocity?) (friction factor). This forces the 
adoption of replacing Kpp by: 
f KpP + በጭ 

The absolute value sign is necessary to make the symmetric p? 
function an anti-symmetric function. Since the KgpPP term is 
primarily due to transverse flow meeting the tremendous flow 
obstruction of the bilge keels, and not for most cases much 
influenced by forward speed, it was decided to let this quantity 
remain independent of ae 

Another form of energy dissipation by rolling or swaying 
motion is found in the formation of surface waves which radiäte 
outward from the ship. This energy dissipation has been esti- 
mated, using potential flow theory, for simple ellipsoidal 
bodies which resemble the underwater form of ships, and has 
been found to depend on ship speed and the frequency of motion 
(assuming harmonic motion). As would be expected, sway gener- 
ates larger waves for an average hull form and ee this 
effect is more pronounced for sway. Thus, it was decided to 
መ Climate sway damping as: 

(Xgy1(o,U) * Yy2(U))v; 

where Yv} is the component associated with wave formation and 
Yvo is the component associated with the hull acting as a lift- 
ing surface in the horizontal plane. Roll damping may be well 


enough approximated by: 
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KpP + KpgPA; 


where Kp is calculable by strip theory, or may be estimated 
om Appendix B as about .05 times critical damping. Kopı can 
be estimated from friction ር tables at an average Reynolds 
number. If the above seems inconsistent, it should be recalled 
that roll damping is much less than the amount of sway damping 
But urs coupled into roll by the equations of motion. Thus, 
the extra care is taken in the estimation of sway damping vice 
BOUT damping. 

TO determine the above coefficients exactly requires that 
extensive model testing be conducted. Since it is envisioned 
that the method proposed by this paper is done at a stage of 
design prior to any such detailed knowledge of the SES. the 
coefficients must be arrived at either by empirical data 
derived from previous similar designs, or else by estimation 
from their theoretical values for similar geometric forms. An 
appendix to this paper gives detailed examples of how to make 
these estimations. 

‚As far as the tank non-linearities are concerned, the 
restriction on the maximum tank angle, n, by the boundaries of 
the tank is the most important of the non-linearities. The 
coupling parämeters Sänd S" remain almost constant for any tank 
geometry as long as centerline symmetry is preserved. However, 
a saturation effect is observed whenever the fluid in the tank 


meaches the tank top because at that point no more fluid transfer 
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athwartships is allowed so that the tank moment reaches a maxi- 


mum. To account for saturation in the equations of motion there 
are two methods. An equivalent linear 'spring constant' may be 
used; to do this change S' such that the same energy storage 

for a given tank angle is achieved for a fictitious unsaturated 
tank as is realized in the real tank. Another method is to 
change the equations of motion at the point that the tank angle 
saturates. After saturation, consider that n = 7 = 0. The 
normal tank dynamics resume again when the fluid can run 'down- 
EUN. This condition is met in the static case for n + g » O. 
In the dynamic case, the dynamic head must be included along 
with the pump pressure, if present. The condition for resumption : 


of normal tank dynamics then is: 


2pAgR2Dv - pAgR2(S"D^ * 2g)f - 2pAgR?gn 


* AgRp(t) » 0. 


The tank damping parameter, By, of course has non-lineari- 
ties associated with it due to the nature of the damping devices 
used in the tank. However, as mentioned earlier, because we 
Enc Dterested in tank design, B+ can be considered a linear 
Coefficient independent of the other variables with the stipu- 
lation that the amount of linear damping must be later designed 
into the tank by some means. It is reasonable to expect that an 
amount of damping in the same order of magnitude as ehealieatrers al 


damping is easily achievable. 
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8. Excitations 

In order to keep the exciting function one-dimensional, 
it was arbitrarily decided to reference all the excitations to 
the waveslope. This was mostly done to be able to compare 
results with previous work. The waveslope used was not the 
surface waveslope, but an average value which the ship senses 
(remembering that surface waves are a local phenomenon which 
decay exponentially with depth). This average value is often 
Fen at the half-draft of the ship, but to be consistent with 
hydrostatics as well as hydrodynamics, I chose to use an aver- 
age value as the value at the center of hydrostatic buoyancy, 
Bach is nearly the half-draft. Also, this places the ex- 
Citations at the center of buoyancy, and in accordance with a 
previous assumption about the non-coupling of sway and roll at 
this location, this neatly separates roll and sway excitations. 

It is now necessary to review linear wave theory for 
simple harmonic ocean waves in deep water. A later assumption 
to be made is that the real ocean is built from a spectrum of 


these simple harmonic waves. 


particle motion | ea Surface ele- 
I S vation 
t 
tv dues Horizontal wave 


a velocity 


Í VU ኮደ. Horizontal wave 
acceleration 

+ e 

AAA ጨ al Waves lope 


Figure 3. Wave Profile 
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Figure 3 is arrived at from the equation of surface ele- 
vation of an harmonic water wave, given in reference 3, which 
corresponds to a surface of constant (atmospheric) pressure. 
Horizontal wave velocity, acceleration, and waveslope are found 


by taking the appropriate partial derivatives. 


[ 


E ea Cos (ks-wt) 


3e/ət + 90° I =-we,„ cos(ks-ut) 


< 
< 
| 


Vw = 9Wiy/9t --u2e4 sin (ks-ut) 


9ፎ/95 = -ke,„ Sin(ks-ut) 


< 
[ 


An important factor to notice from the above math is that 


at any given point in the wave y and v, are in phase with each 


ሃሃ 

other but are out of phase with the wave horizontal velocity vy . 
For simplicity, the wave variables may be expressed in 

complex form, with the implication that the quantity indicated 

is the real part of the complex expression. 


e — ea exp(1i(ks-ot)) 


NN ~ E 0) Vw = iw*e; ሣ = ike 


Vw/ 후 -w/ik = ig/w 
VWw/y = we/k = g 
Mae last two relationships fix vw, and v, as functions of y. 


Assuming that the excitations occur at the center of 


buoyancy, and using the same prescription as before: 


Ko (t) — K'yy + K'y vw * IEA 
ነር!) >"ሮ፡ + Y y Vw + Y' yu VW 
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There is an unrealistic assumption in the above derivation, 
however, in that the horizontal wave motions are not impinging 
on a static ship. They are actually influencing a moving ship, 
which suggests that for the equation to be valid, v, and ቀ 


' and = - v, respectively. 


Mist be replaced by w - v 
By making the assumption that the hydrodynamic coefficients 
of the water moving past the ship are the same as for the ship 
moving past the water, and by adopting the classical assumption 
of Froucethat Kj, = Kg, all of the coefficients are determined. 
Notice that Y || vanishes due to a previous assumption, but that 
K'y,, and K' cannot vanish unless the moment arm, zp , between 
the center of buoyancy and the center of gravity (the reference 
for the x-y-z axes) vanishes. With all of the above assumptions 


and conditions included, the excitations may now be found. 


—- V 7 212 


VES zh D 


ia? = 84” መፕ yu (VW ነጀ! Zhp) + zhYv (vw-v F zhp) 
YO(E) = YXv(vw-v t zhp) - YV (ww - v + zhb) 
All of the terms in the excitations, except for those in- 

volving U,Vw(V) and DE may be moved to the left-hand side 
of the equations of motion when introducing this expression for 
the excitations into the equations of motion. When combined 
with like terms, this doubles the value of some of the hydro- 
dynamic coefficients, and introduces coupling terms into some 
others. Also, the excitation expression added two new parameters 


instead of the one desired--the waveslope, and the phaSing 
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difference between waveslope and horizontal wave motions. 


E Equations of Motion 


The results of the previous sections are summarized to 


completely form the equations of motion for an anti-rolling 


Bank. 
Roll) [ (Kp-Ix) D2 ID “10 - zh[YvD + Yv] (w,U) + Yv4(U)]v 
- pAoR? (S"D® + 2g)n + Kagel - Kyß? = -Kolt) 
(sway) zh[YvD^ * (Yvj(o,U) + Yvo(U))DI¢ + [(¥v-m)D 
+ Yvy (w,U) + Yv2(U)]v * 2eAgR^D^n - -Yo (t) 
(tank) RAD 20) MEE A RD 
-pAoR? (S'D2 + (Bt/pAoRZ2)D + 2g)n = 0; or- AoRp(t) 
E uwreetor) (TpD + l)p(t) = queer (G5 D) Cute ( geet om DV 
MS PEN OE se 
Ecuacrons) Kolt) = -Kgy + zh(Yv] + Yv2)(vw - v + zhp) 


+ ZhYv (Vo ፦ v + zhp) 


Yo(t) 9 -(Yv4 * Yv3) (vw - v + zhp) - Yv(vy - v + 
zhp) 
where: n = n = Oand|n| = nmax, whenever Inlexceeds nmax. ; 


the condition for returning to n and n free is 
whenever: 


2 
2pAoR^py - 0 AoR? [ s" p2 + 29]# - 2pgAoR2n * AoRp(t) » 0; 


Wiehe ele whimaelab condition n = + nmax and o 
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These equations may be linearized by dropping the terms 
involving 2 ao, and the restriction on tank angle and 
pump pressure. A more convenient way to express these equations 
is to non-dimensionalize them; the results in non-dimensional 
form represent a family of solutions which help to give an in- 


sight of the general design problem. 


10. Non-Dimensional Equations of Motion 

Because further analysis in this paper uses mostly the 
linearized version of the equations, only the linear terms are 
given below. A paragraph at the end of this section gives the 
non-linear terms as well. 


Dimensionless variables were chosen as follows: 


US SE 
ED “/ወሠርጅ 
ከ. p(t)2pgR 

 ዜየከ ww in radians 
[) = dl )/dt 
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The linearized equations, with complex notation for the 


K Citations, are: 
(roll) [(1 + 1p? y1) D? + (Kg + ApéKy)D + 1]2 
E 2A, ly D + Ky] v* ate At (D2/ ust? + l)n 
= [(1 - Apdy) ~ iàìpàìo] v(t) 
(sway) -2/b[yjD? + KyDJ@ + (Y5D + 2Ky) v* 
_ ED? / 가, = (ày + 120) wc) 
(tank) AE(D4/ug42 * 1l) - AtD/uye2 v* 
* At(D?/uc2 + KtD/ut + 1)n = 0; or At p* (rt) 
(effector) 209R(wsTp + 1)p* = (g] T ws92D)2 + (93 + wsg4D)n 
Coupling constants Damping constants 
At --20gAoR^/Kg Ks = Kp/ VKg(Kp - Ix) 
Ab = Ky = B2YyV Kg (Kp - Ix) 
E = Kt = Br/pAoR2 V2g S' 
Eo -—(Yv] t Yv2)gB/uoKg 
Enertra constants Frequency constants 
Al = B?Yv/(Kp - I4) ut = wt/ws 
Ao = B2 (2Yy 7 m) / (Kp zul) Ust = ሀ 5ቲ/ ሠ 5 
Uyt = W9/B/ws 


Spring constant (ship) 


Kg = -mgGM 
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In order to preserve the non-linear equations, the follow- 
ing terms need be added to the indicated equation, and also the 


restriction on tank angle and pump presSure applies. 


0011) Kappl + 303 


where: Kg lus? Kop) /ቾጋ 


፪2 = 1/2 , or l/ég? 


There 1S one great idiosyncracy of these equations which 
would pass unnoticed were it not for the presence of the wo term 
Ehe wave velocity - sway coupling constant, 45 . Here, wo 
refers to the absolute frequency of the waves relative to the 
Sartn. The rest of the dynamics will occur at a steady-rate 
frequency which is the encounter frequency that the ship senses 
in the waves. This encounter frequency is a function of both 
ship speed and relative heading angle, and is further examined 
mae the appendix on excitation spectra. The peculiarity is that 
Ehe ship is sensitive to both relative and absolute excitation 
frequencies, which is a direct result of the choice of sway 


excitations. 
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NER. Excitation Spectrum 


l. General Discussion 

The preceeding chapter completely defined the anti- 
rolling tank problem with the exception of the single inde- 
pendent variable y (time). This chapter is dedicated towards 
obtaining a realistic expression for Y. Ocean waves are 
essentially a random process. To account for the statistics 
of the wave process it was necessary to adopt an accepted 
empirical spectral (implying frequency) distribution function 
which had the disadvantage of being an elevation distribution 
instead of a slope distribution. The main part of the text 
of the chapter deals with making corrections and additions 
>> His distribution function. 

It is first changed into a waveslope distribution by use 
ec linear gravity wave theory. Correction factors are derived 
for the amount of submergence of the ship's hull, for the 
relative heading of the ship with respect to wave progression, 
for the wave encounter frequency shift--akin to the familiar 
Eus. shift--and for the shape of the ship's hull. Because 
er the highly speculative nature of some of the assumptions, 
particular attention is als to their expected range of 
validity. All of the factors are then combined to form a 
complicated expression for the excitations which can be 
inserted into the equations of motion to complete the model 


of the ship-tank system in ocean waves. 
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This strategy was adopted because it lumped as many of 
the parameters which are independent of the anti-rolling tank 
design aS possible into a single expression for excitations. 
mies could then be used as an input forcing function. The 
resultant form of y- (time) had the peculiar but advantageous 
properties of being "steady state", i.e. a combination of 
regular harmonic waves, and linear while only requiring that 
the various auxiliary component parameters be analytic. Thus 
the final derived expression for the excitation spectrum y 
included most of the ee independent parameters as well as 


the statistical distribution of waveslope. 


2. Waveslope Distribution 

In order to describe the wave excitations of the sea 
surface it was decided to use the empirically derived Pierson- 
Moskowitz wave amplitutde spectrum. This is a standard 
analytical formulation which seems to adequately describe the 
real life wave amplitude spectrum for fully developed seas. 
The statistical distribution of wave amplitude over the band 
of wave frequency (0 to ~) is assumed to depend on just one 
parameter, significant waveheight, ma Significant wave- 
height is commonly used as a measure of sea state. To be 
more specific, EUN wave amplitude squared per unit fre- 


quency bandwidth, is a function of frequency only for any 


given sea state. For a more complete description of this 
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spectral density function, its implications and proven uses, 


፡ ቲ ር | ፦፤ references 3 and 11. 


c am = 0.0081 com ee) 


1/3, 2 


W 
| 


0.0324 492/ (8 


32 ri ee 


Ko) 
I 


Using linear wave theory, it was previously shown that: 
EINEN s] = ke 
2 
where: k =w /g 


Making the appropriate substitutions, the spectral den- 
sity function for waveslope is found. 


4 
U^ /6u N < 8, , (0) _ ug e le 8/0 ) 


The mean square value of the waveslope spectrum is found 
by integrating the spectral density function over the entire 
frequency range. 

22 co 


_ 2 
V = K S e g (w) dw 
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The root mean square value (RMS) of wave slope is found 


/ 


by በከር. = Ü. By the statistics of a random process with an 
assumed Rayleigh distribution, the average of the 1/N highest 


values of y are given in reference 1l: 


yMN . የ(1/8) ፡ ህ 


EMS 


1/N £ (1/N) 
1 NEL 
1/3 2.00 
1/10 ን 
1/100 3.34 
2208. 3.29 


It should be pointed out that the above numbers refer to 
Ehe P ገ0ርሬ angle measured from the horizontal, and not to 
the "crest to trough" condition usually associated with 
waveheights. 

Following the convention that "Significant" means the 
average of the 1/3 highest values of the process, the signi- 


ficant waveheight is: 
S 
Y l ር 


This represents an acceptable peak value, for engineer- 


ing purposes, of the random waveslope process. Notice that 
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the mean square value has a different relationship in this 
random process than for the deterministic case in which peak 


value is 72 times the RMS value. 


B. Depth Correction 

The surface elevation of the sea is a constant pressure 
surface. This pressure decays exponentially with depth until 
-at a certain depth, the surface waves are not felt. By linear 


. : ! -w*zo/g 
theory, this exponential decay factor 1s e ; 


where Zo 
is the depth from mean surface elevation. 

Since the center of buoyancy of the ship is the point 
through which wave excitations are induced in the static case, 
and for small amplitude excitations (i.e. nearly static case), 
a reasonable approach towards excitations is to apply them 
at the center of buoyancy. The only alternative to this is 
to carry out the integration of the pressure distribution 
over the entire immersed area of the hull. By the above 


reasoning, the spectral density function which the ship feels 


may be approximated by: 


-2tbw*/g y 
e E (0) ; 


where ty is the draft to the center of buoyancy. 


Bees COrrection factor for depth is shown in figure 4. 


Merce that for super ships, with drafts of 80+ ft and t, of 


b 


about 40 ft, waves are practically not felt in roll 
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ር  ተ1ነ). For larger ships, say 20,000+ tons, all wave 
ከ ዎክ less than about 50-100 ft. are unfelt, and are non- 


interactive with the ship. 





20,000' —200' 50.6' 22: 50 
Figure 4. Depth Correction 


4. Encounter Frequency 

Since the ship forward speed and wave advance are both 
vector quantities, relative encounter frequency of the ship 
with the waves changes with orientation. Denote heading 
angle of the ship into the waves by B (B is 180? in head 


Seas). The applicable relationship is: 


wo = w (1 = WU cos B/g) 
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5. Transverse Waveslope 

The equations of motion were developed for roll and sway 
coordinates which have entirely a transverse sense with 
respect to the ship. It can be shown that for waves progres- 
Sing at some other direction than transversely to the ship 


the effective waveslope in the transverse direction is: 
UR) 61006); B is 90° for beam Seas. 


Region Of Validity of Excitation Spectrum 
Combining all the corrections to the original waveslope 
spectral density function, an expanded function is formed. 


_ 2 
y?/ou(w, t,, B) - sin?(B) e ?*b" /9 g.. (u) 


EC 

The expanded form of spectral density is plotted on 
Figure 6 for beam seas (sin(8) = 1) for the normal range of 
expected sea states. A glance at the figures shows that for 
all ships with a draft of 5 feet or greater, the desired 


mange of integration of S ር ንጋ. 5ጤ.፡፡፡ 4 1|!!:።>ከክ 


EC 


s The limits on w impose the corresponding wavelength to 
fall within 20 < A < 20,000 feet. 
To determine the range of Ë possible, a visualization of 


the ship in long crested waves is helpful. 


내 
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i 
wee =  — - ዴዴ. 


ie a cum 
y U 
tg aa ሪ ዜ dl 


| L|cos (A) 
ae ae 


Bagure 5. Ship in Long Crested Waves. 
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Significant Wave Height = 30 ft. 


.01 Beam Seas 
yl/3*= 30° 
y 
tb=0 
5 
10 
0 | 
1 2 a. TS 
01 Significant Wave Height = 15 ft. 


Beam Seas 


y1/3 = 15' 


+ 
Ku 





zo Significant Wave Height = 5 ft. 
Beam Seas 


yl /3=5 | 


I 





Figure 6. Excitation Spectral Density in Long-Crested Waves. 
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It seems reasonable to assume that the excitation at the 
center of buyoyancy is equivalent to that for the whole ship 
as long as the value of L|cos(ß)| is less than a half- 
wavelength. Another consideration is that the roll excita- 
Buono) the more significant of the excitations, is propor- 


Euonal to Kg. A look at the terms in Ky is enlightening. 


Se = pg(Immersed Volume) GM 
E — (a constant) GM 
BM = (2nd moment of waterplane area about x-axis) 


(Immersed Volume) 


N 


GM BM + KB - KG 

The above implies that x 1S proportional to the 2nd 
moment of the waterplane area. It can be shown that for a 
rectangle, the 2nd moment is Therefore, it can be 
seen that Sg iS proportional ae DE 

Since the beam is not constant along the ship's length, 
End, sthe wider portions of the ship will contribute to roll 
excitations. It is easy to calculate that any section with 
less than 60$ of the maximum beam will have a negligible con- 
tribution. Therefore, the ship can be considered to have an 
effective length Le which should be used in the formula 


L|cos(8)| to replace the waterline length L. An empirical 


formula, which is a conservative estimate based upon the 60% 





39 


max. beam criterion, was derived graphically using a Taylor 


Series hullform from reference 9. 


+ 0.943 Cp 





Figure 7. Effective Shiplength. 


The region of validity for heading angle is therefore given 


by: 


66 E IA 


DUC ss- A = 27 


Applying this standard to graphical display is done in 


Figure 8 below. 
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In-valid 
Le 600' 


0.5 1.0 1.30% 


Figure 8. Region of Validity for Heading Angle. 


To summarize, it is obvious that the range of parameter 
values for w and B vary significantly with ship dimensions, 
ship geometry, and sea state. It is important in each case 


to investigate whether the extended equations are still valid. 


7. Heading Angle Correction 

From the results of Section 4, it would seem desirable 
to extend the range of heading angles further by some approx- 
imate means. An easy way to do this is with the assumption 
that the static righting moment is uniformly distributed 
along the effective length of the ship. Then, the wave exci- 
tation vanishes anytime that Le|cos (8) |/ has an i (Cru ji 
Of 1 or greater. It can also be surmised that for non- 
integer ae ate cos (EN gredtermwbhan i» the part of 
the ship's effective length between successive crests 
receives a net zero excitation due to cancellation effects. 


Tt is appropriate to include this effect as a correction 
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factor to the excitation spectrum. An efficiency parameter, 
€, which is a function of Le, B, and A(w) is formed by put- 


ting the above arguments into mathematical form. 


et ከ ር À 


X is between e[Le, B, X00) 
0.0 0.5 1.0 
0.5 መሙ 2(1 - X) 
I0 1.2 2/3(X- 1) 
T: 5g» DES D - 2) 
2 በመ... ጋሪ ን በደ. "=" 2) 
2.5 3.0 2/5(3 - X) 
1 5 2/7(X - 3) 
60504 . 0 2/7(4 - X) 
420 425 «0 
€ 
€ N 
나 
1/3 
E 
ge vA 1/9 
0 1 2 3 4 X 5 


Figure 9. Heading Angle Correction. 


4 
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The excitation spectral density function now is: 


2 
V / 8 = El 518" (6) eee ee (w) 


Se 
Since an actual ship will probably not have a constant 


beam along its effective length, K, will not be uniformly 


ቀ 
distributed. The region of net zero excitation will not be 
quite equal to one wavelength, but will change with the 
relative orientation of the passing wave. Qualitatively, 
enis will tend to M dice higher harmonics into the 
encounter frequency of the roll-inducing wave energy than 
that predicted by the above model. However, in general it is 
safe to assume that the cancellation effect does occur, even 
though the higher harmonics are produced. 

Ku L her Consideration is that the final product of 
interest, roll angle, will be determined from integrating the 
spectral density function, and taking the square root of the 
result. This set of operations makes the roll Me very 
insensitive to the exact form of the heading angle correction. 

Based on the above comments, it seemed consistent with a 


conservative approach to take the correction factor as the 


line connecting the peaks of the function shown in Figure 9. 


1/2X, otherwise. 


ሰ 
| 
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8. Spreading Function 

A spreading function was used to convert uni-directional 
long-crested seas (implied by Spp! into a more realistic 
short crested sea surface. Following reference 3, the 
e anng function f(u) is a symetric function centered on 


the primary direction of the seas. 


f(u) = 2/T cos? (u), -7/2 < ህ < 1/2 


/"/5 ቶ (ህ) 6 


(Jdu = S 
7 56 


W 
gg No) 
The effect of this function is to smooth the long 
crested seas prediction for roll response somewhat. It implies 
that there is always a beam component of wave motion influenc- 
ing the ship even when the ship is headed into or away from the 


predominant direction of the waves. 


9. Sea State Data 
The below listed values were taken from references 8, 11 and 
12, and are useful for design purposes. The peak frequency 


of the surface wave elevation spectrum S,, is the መ 


።: 
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17: Probability of Occur- 





0 DRM Sea State ence in North Atlantic 
02417 30 8 0.01 
02455 25 7+ 0.02 
0.509 20 7 0.03 
0.589 LS 6+ 0.10 
02720 10 5 0.35 
1 . 0 2 05 3 0.45 


ር ትጋ <05 O 0.04 
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fee Solution for the Equations of Motion 


fo Introduction 

Given that the coefficients are all known, there are at 
least three practicable methods of solution available. These 
are analog, numerical integration, and complex-algebraic 
manipulation of the linearized equations. Each of these 
methods has its own advantages and drawbacks compared to the 
others. The latter method was the one finally adopted because 
it was the most flexible for design purposes and yields the 
most useful information for effort put into uSing it. How- 
Ber, Partner work with the other methods may be warranted 
depending upon the particular anti-rolling tank problem being 
investigated. The details of the complex algebraic method 
are given in this chapter. A brief overview of the results 
of the investigations into the other two methods is given in 


Appendix D. 


E Complex Algebraic Solution 

The use of complex algebra along with the spectral den- 
Sity technique allows for a closed form solution of the 
linearized equations of motion. The basic assumption involved 
in this solution bears a close resemblance to what may be 
imagined as the real situation of the response of a ship in 
waves. It should be recalled that gravity waves have an 


advance velocity which is dependent upon frequency, so that 
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the phasing between ocean waves is essentially random. The 
waves of different frequency (or wavelength) travel along at 
different velocities. If it is assumed that the ocean has a 
constant spectral density with time, then the ship's response 
to each of these frequency components has reached a steady 
state value. Random phasing of these components accounts for 
the Rayleigh type distribution of peak values of response, 
and the principal of linear superposition neatly ties the 
whole system of equations and density functions together; the 
ship's response is the sum of the steady state responses to 
each of the frequency components. 

A mention should be made of the implications of the term 
"linearity". The basic mathematical operation for each com- 
ponent case of frequency excitation and response is linear 
because the coefficients have been determined. However, the 
E 41685 for the so-called linear coefficients which 
involve Qr Wor X, phasing, ship dimensions and geometry, 
tank saturation, heading, sea state, etc. have no such 
restrictions. They are required only to be analytical. 

The standard method for finding the steady state solu- 
tion to a set of ordinary linear differential equations is to 
make a complex transformation, using a complex algebraic 
Operator for differentiation, solve the complex equations 
algebraically, and retrieve the information of response mag- 


nitude and phase angle from the complex expression. 
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Let Y = ር ከ Then the steady state solution must be 


Àt 





, a iut u k = wk 
ethe form: 0 vac 86) ካዕ es V KS , where 
* 
i. $o: ng: Vg are complex constants. The operator 
D = d( )/dt may be replaced by D'= iu, p? = E The equa- 
tions of motion now become for a passive tank: 
.1!1 (1 ነከ8ከዜ ቁ... ብ) ር 01711111 ue LAA 
: ! እ ብ i 
er OT .(ህ065(.273(114) ue un Vo 
qq MEN 3209 C33 (111) no 0 
Qe l - Ay A + IA AG 
C (iu) vo = = + IA Y 
no 0 
N (iu) 











do 
DS 


determinant [Cí(iu)] 


The numerator, N(iu), is found by replacing the first 
momomm Of the C matrix with the coefficient matrix for Vo - 
Notice that the C matrix terms will not be symmetric because 
of the choice of v* as a variable. If y* were chosen instead 
of v*, the symmetric condition for an energy conservative 
System would have been met. 

Since p = 110 by linear theory, the roll velocity 


response may be obtained by uel. and likewise the roll 
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acceleration is gel. Tank angle, velocity and accelera- 
tion are similarly determined. 

The final step in obtaining roll response is to multiply 
the transfer function ġ?/y? (or àó ?/y?, or $9 ?/y?) by the 
EUrtatjijon spectral density function 1*/6w and integrate the 
result over the entire frequency band. The spreading func- 
tion f(u) may be inserted in this step depending upon whether 


long-crested or short-crested sea result are desired. 


mee c - 
De sa Y / (Ol ou cO (long-crested seas) 
1⁄2 VE on an 
ome = 4 FC) /S (6*/0*) 07 /du dwdu (short- 
n/a 0 crested seas) 
Maine 


ris 
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V. Anti-Rolling Tank Design 


I Introduction 

The complex algebraic form of solution of the non- 
dimensional equations of motion was written in computer 
language to allow the numerical work to be done by a digital 
computer. Details of this computer program are given in 
Appendix C. Essentially there are two separate programs-- 
MAIN 1 which maps as versus M and Ky for one case of head- 
ing angle in long-crested seas and MAIN 2 which calculates 
short-crested sea results for any given set of input para- 
Meters. In the interest of numerical efficiency most of the 
iterative steps of the design process can be done by MAIN 1l, 
looking at the largest component of roll in long-crested seas 
which makes up the center frequency of short-crested seas 
results. MAIN 2 can be used later in the final stages to 
give short-crested sea roll angle predictions based upon the 
results of using MAIN 1. | 

The design method basically entails "bouncing" back and 
forth between MAIN 1 and MAIN 2 looking for optimum para- 
meters of the non-dimensional equations. The steps are man- 
ually done, meaning that each step represents a separately 
submitted computer run. About 2 to 4 computer runs are 
required for a final design for each different anti-rolling 


tank. 





50 


The passive tank design method in section 2 suggests how 
to specify various parameters of the equations and presents a 
logical sequence to arrive at an approximate tank size. 

There is also a suggested iterative sequence to be used as a 
design aid. One of the fortunate results of the subsequent 
chapter on design examples is the collection of equations 
shown on Figure 12. It is recommended that this worksheet be 
kept as a tally sheet and guide when performing the iterative 
steps. 

The active tank design method given in section 3 is 
mostly a discussion of possible changes to the passive tank 
design which are compatible with the computer program. Mak- 
ing these changes will result in a prediction of the perfor- 


mance of a feasible activated anti-rolling tank. 


2. Passive Tanks 

The response of the ship-paSsive tank system is modeled 
by the non-dimensional equations developed previously and 
setting F (t) = 0. 

Tank location may be specified or a small number of dif- 
ferent locations may be investigated. The reason for limit- 
ing tank location to the absolute minimum is entirely practi- 
cal. The vertical space of the ship is divided by structure 
into a number of decks and for access and strength considera- 
muons the cross duct of the U tube tank should run across (or 


under) a deck. Also there may only be a single or a few 
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Bons within the ship large enough to accept the tanks. As 
a consequence of the non-dimensional equations, there is a 
curious coupling between tank vertical placement, size, and 
the coefficients of the equations. Since the ship's vertical 
center of gravity is shifted with any vertical shift of the 
tank, the ship's GM is reduced whenever the tank is moved 
higher in the ship. Since E is proportional to GM, and We 
is proportional the square root of Kae and because all of the 
non-dimensional coefficients depend on ES or wo, each verti- 
cal tank location must be treated as a separate problem. The 
design method which follows is for each separate case of ver- 
tical tank location. 

Once tank location is specified, the parameters which 
NM Lu ሺ. We know 
that Ky has a minimum value of about -0.05, and that M will 


have a value close to 1.0 because the vibration absorber 


ከክ ክክ. ፡]] free to be chosen are: A 


effect of the tank is maximized by tuning the tank to 


uU, = w/w = 1.0. These four dimensionless parameters are 


e S 


dependent upon the real ship and tank parameters as follows: 


È ጋ 

해 = A, IR, (CM) , PA R ‚m] 
Um = ይ. R,S' (A,/Ai ,L. /R)] 
L Ho, [ሠሊ (ጆ ) R,S"(Z /R,L /R)] 


_ 2 
Kt = K [pA R , geometry of tank, A/A) 
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Berne ship Changes its loading condition during 
expected operations such that its displacement is appreciably 
changed, then all four parameters have a A dependency, also. 

Another important parameter for design is max. Tank 
capacity y max is defined as the maximum heel angle that the 
tank can impart to the ship by maximum movement of the tank 
fluid off to one side. This also represents the maximum 
waveslope angle which the tank is capable of stabilizing. To 
form an opinion of a required yy max, use must be made of 
empirical formulas. Hagen mentions the use of: 


D a 0.36/Log,g (A) radians 


A formula suggested by center of buoyancy correction to the 
spectral density function which takes the equivalent of a 6° 
surface waveslope at the center of buoyancy for the ship's 


natural resonant frequency is: 


E 2 
y = | ) e (tpus*/9) cae 


ax 
The above formulas give an idea of the approximate capacity 
which a tank must have in aa to be effective, and serve as 
a feasibility check before proceeding further with a particu- 
lar design. 

The parameter A, = (~2pgA_R°/-mgGM) has certain char- 


teristics which allow it to be specified by the designer. 
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Rearranging the formula to read A, = [2A /(016018000 volume 
of ship)]/(GM), notice that the numerator of this ratio is 
the standard free surface correction for a virtual rise in 
the ship's center of gravity due to the fluid being unconfined 
Ma the n direction. Therefore, Ar is the percentage loss in 
emne stability of the ship due to the tank fluid free sur- 
face. There will be a maximum allowable hy which ensures 
safe operation of the ship under all loading conditions, 
including damaged conditions. Since a larger A 1S associ- 
ated with more tank stabilizing effectiveness (Pax ¿ER y 
pne choice of Ar is fixed at the maximum value which allows 
safe operation of the ship. 

An important conclusion to reach from the above observa- 
tions is that the ship hull form should be designed with the 
notion that an anti-rolling tank is to be installed. This 
means that space, weight, and static stability will be such 
that the ship will be comfortable with a tank of the capacity 
indicated by the empirical formulas. The alternative choice 
Sie trying to fit an anti-rolling tank into an existing ship 
design which may not have all the necessary reserve space, 
weight, or static stability most likely will result in an 
ineffective stabilizing system. 

Now that AL is known, and ER has been determined 
approximately, an assumption for a Hr value is all that is 
required to calculate tank dimensions (Ayı R, Lo: Z.) the 


tank parameter Hoy (LQ/R,2Z,/7R) , and approximate tank weight 





54 


[2 pA R (Lo/r 4 A/A )]. A first guess at u, is made by using 


는 
Chadwick & Klotter's equation, based upon the linear equations, 
neglecting sway, from reference 4. 


B h 2 od A 2 
m [1 Ay UU 1/ዞ..+ ) ] 


It is now possible to calculate all the tank dimensions 
and all the coefficients of the equations of motion with suf- 
ficient accuracy to consider them constant, except for duct 
area A4 and tank parameters M and KL: 

The remainder of the design problem may be stated as 
follows: Find the best combination of M and K. which mini- 
mizes the rolling motion predicted by the non-dimensional 
equations of motion. This two parameter search problem may 
be solved by noticing that Hr must be about 1.0, and that two 
it accuracy is the best to be expected from the equations 
of motion. Therefore, instead of a sophisticated search 
technique it should suffice to map roll angle versus Hr and 
ሺ. This may be.further simplified by sub-optimizing each 
Wy value for Kye This approach is justified by gentle cur- 
vature of the function ON in the region of interest, 
Euound in all sample problems tried. The remainder of this 
section will investigate numerical methods to deal with this 
search problem. 


The basic strategy for solving this problem is to first 


fix upon some criterion of roll performance. Based upon the sea 
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State data already given, the author decided to minimize roll 
angle at the worst heading at operational speed in sea state 
SMS feet significant wave height) was a reasonable criter- 
Meme Other criteria could include roll velocity or roll 
acceleration for other environmental conditions. 

In order to compute E for given conditions, a long 
sequence of operations is ee 1s twice the 
square root of the integral of E times ay cae This 
sequence makes vivid the reason why the complex algebraic 
method of finding $/U is superior to numerical integration. 
The latter would require another integration of many steps, a 
E expensive approach since the results must in turn be 
iterated many times for a good design. 

The above discussion assumed that long crested seas were 
involved. For short crested seas another integration, this 
time over every possible heading angle, must be performed for 
the spreading function. For this reason, long crested seas 
are assumed for mapping o K, (best)]. Short crested 
EEuscare only looked at later as an output. 

In mapping el qu, K, (best)], it was found that the 
function displays very gradual curvature as was to be 
Expected for an integral function. Therefore, only a few 
points are investigated during each computer run, because 
meer polation is predictable. A second order Legendre 


polynomial curve is used in an interpolation scheme to find 


a local minimum or maximum. Manual inspection of the computer 
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results will show whether further mapping is necessary. 
Warning: the interpolation function likes to extrapolate to 
find minimums or maximums equally well, but such results are 
invalid by the statement of the problem. 

Once iu K, (best) ] has been mapped and values of 
M and K. for a minimum roll angle have been decided, short- 
crested sea results may be obtained. The problem then shifts 
to chosing which of the alternative tank locations is the 
best. Economic as well as performance factors will enter 
into the decision of this question. For example, a tank 
which is in a cargo hold, presumably valuable real estate, 
would probably not be chosen over one which was in a less 
costly part of the ship, such as atop an outside deck, were 
the two tanks to perform nearly as well. 

A detailed listing of the author's steps in passive tank 
design are reported as an aid in improving upon them. A 
later section will include a completely worked example of the 


method. 


I. Problem setup 
a) Assume a range of a, small number of vertical tank 
locations. 
b) Calculate Wa requirement. 


c) Specify a እህ 


d) Assume an initial Hy value from Chadwick and Klotter's 


equation. 





UT. 


UII. 


EV. 


e) 


£) 


g) 


a) 


= 


Determine Lo, Aor R; 247 and approximate Ay - 
Design the tank except for A, - 


Calculate the hydrodynamic coefficients and other 


Neomputer Inputs: 


Brest Iteration 


Run MAIN 2 program; use the long crested wave results 
to estimate the worst speed and heading angle com- 


assumed, K, = -0.8 


bination. Use y a 


t 


b) Map $3, K, (best) J using MAIN 1. Use Hy assumed, 


C) 


K initial = -0.4, ለ9. 


Interpolate to find ad 


0o AR, = -1.0. 
73 


minimum. Check computer 


Beeuleszeosensure that a true minimum occurred. 


Second Iteration 


a) 


b) 


Run MAIN 2 again to see if the worst speed and head- 
ing angle combination has changed. Also, roll per- 
formance may have degenerated at other headings and 
speeds significantly. 

Continue iterations if change in speed or heading 


angle warrant it. 


Alternative Tank Choice 


a) Based upon the operator's or owner's criteria, select 


the best tank location. 
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Ve Verification 
As a final step, short crested sea results, MAIN 2, 
should be run for the entire expected range of speed and 
sea conditions. This is to make sure that the solution 
didn't fix too shortsightedly upon some assumed worst 
condition and lessen it without overlooking its negative 
effects on some peculiar or unforseen combination. 
Another verification which should be performed during 
this and all other steps as well is to check the range of 
parameters, the range of the solution variables to ensure 
that their coefficients remain linear (or else adjust the 
coefficients to their equivalent linear values) and that 


the M and K, values being predicted are realizeable. 


t 
E Acti)vitated Anti-Rolling Tanks 

An activated anti-rolling tank is one which has an input 
power source in addition to ship induced excitations. A 
pressure source is used as the type of power input in this 
discussion, as opposed to a flow source. It is desired to 
regulate the pressure source in such a way as to reduce the 
melting motion of the ship. 

The choice of a pressure source instead of a flow source 
was not entirely arbitrary. A flow source would have to be 
some sort of a positive displacement pump capable of handling 
large volume rates of fluid transfer, approximately 500 tons 


per minute as an order of magnitude. The required pump size 
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for this volumetric capacity would make it economically 
unfeasible. Another requirement for the power source that is 
imposed by economics and amply illustrated in other papers, 
Em reference 18, is that it operate at a nearly steady 
power level. The peak power output for a pressure source, as 
found by standard control theory methods using linear feed- 
back, is too high to be practical. There is too much wasted 
power capacity which is only used at peak power demand.  How- 
ever, the pressure pump must be built large enough to be able 
to supply the predicted peak demand if the linear equations 
are to be valid. Once further economic factor to be con- 
sidered is the competition, the alternative methods of roll 
damping. With regards to these, the gyroscopically controlled 
anti-rolling fin installation is the most serious candidate. 
In practice, a standard, well designed anti-rolling fin 
installation has been found to consume about one or two per- 
cent of the installed hoursepower of a ship, including pro- 
pulsion, due to the added effects of the fin area increasing 
the ship's total hull resistance to passage through the water 
and the power needed by the fin positioning motors. Thus, 
economics and practical engineering considerations have dic- 
mated a rather narrow — possible control effectors for 
anti-rolling tanks--a pressure type source of power of nearly 
constant output which has a power rating of about one percent 


of the ship's installed power. 
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A further consideration of the power source is reversi- 
bility. A level power output does not imply that the direc- 
tion of pressure or flow will be constant. The pressure 
supply will obviously have to reverse itself as the ship rolls 
back and forth. This aspect of the pressure source, quanti- 
fied by a characteristic time lag to change direction, will 
have to be reasonably estimated and included in the analysis, 
and must also be reasonable in the sense that a smaller time 
lag probably represents a more expensive piece of equipment. 

It is helpful to first look at the orders of magnitude 
of the parameters needed to describe a control effector. 
Consider a ship of 20,000 tons displacement with a speed of 
20 knots. Assuming an overall lift-to-drag ratio of 150, the 
installed horsepower would be about 18,400 Hp. This would 
indicate from the preceeding arguments that an economically 
feasible anti-rolling tank can be no greater than 184 Hp. 
Compare this to the energy flow of the ship in its roll 
Seeraanate, assuming a GM of 6", a 5° roll, and a roll period 


of 12 seconds. 


0, 000 eons (2220809 4c60n) (6') (sin. 5?) (1/2) - 


(1/12 sec) (Hpsec/550 ft-lb) = 1775 Hp 


This simple estimate is enough to show that we cannot 
pect the control effector to be able to counter the effects 


Of wave excitations. However, it may be used to counter the 
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resonance phenomenon associated with rolling frequencies near 
the natural resonance peaks of the ship-tank system. 

Since the passive anti-rollsng tank as capable of devel- 
oping nearly as much energy flow as the ship, until the 
saturation angle of the tank is reached, it seems justifiable 
to further simplify the control effector equation to 


Js PEE pt) = an 


Po denotes the steady pump pressure developed across the tank 
Auct. and T, is the reversing lag time for the pump. The 
equation P. = £(6,6,n,n) implies now that the sign for Po is 
chosen according to some criterion contained in the state 
variables $,0,n, and n. Consider the same ship with an 
anti-rolling tank of 200 tons weight of fluid. The average 
fluid volume flow would be: 


807(200 tons) (35 Pe E l O) (1⁄127Sec)m= 142 Er 


p 
For an average 184 Hp, the mean pressure is 184 Hp (550 ft-1b/ 
Hp-sec) (1/142 ft? /sec) = 245 Le E = 1.7 psi. This pressure 
changes the fluid level between the tanks by 3.8 ft. Assum- 
ing a 50 percent pump efficiency, the pressure differential 
would be 0.85 psi with a fluid level difference of 1.9 ft. 
This corresponds to a tank angle of about two degrees. The 


Static ship roll angle for still water and constant pump 
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pressure would be À. times two degrees, amounting to only a 
fraction of a degree. Thus we can now state that the addi- 
tion of such a control effector would have a qualitatively 
negligible effect on the equations of motion, except near the 
resonant frequencies. This indicates that some form of 
active damping is the goal of the control effector switching 
eriterion. 

The problem of adding a control effector now becomes one 
of how to add the active damping to the linearized equations 
of motion. The tank equation is the only one of the set 
Kun contains p(t) for an active tank. Substituting Po -— 
915 + 935 + g4n + gn into the tank equation, neglecting for 
the moment N (consider i small), and combining terms in the 
linearized, non-dimensional equations of motion, the struc- 


ture of the active tank equations is revealed. 








ü 9 g 2 
2 2 s72 l D 
a MM st o 26፡8. E 2 E ~ 
My t 
K ü g g 
2 2 E s74 ws”73 _ 
+ A, [D ^ + m 2 UD t (1 - 2698)" 0 


We would expect that the effects of 91 and ፪3 to be 
relatively small even if we took them at their maximum pos- 
sible values (to develop 184 Hp) because of the small static 
values of tank and roll angle the pump is able to develop. 


Hence, making 9, = 93 = 0 is approximately correct for any 
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Seenomical tank. The choice of 2 and J4 to use as a switch- 
ing criterion is obvious. The tank already has linear damp- 
E. which may be varied as required by adding valves or 
baffles. Negative tank damping is undesirable because it 
would increase roll response at resonance frequencies. The 
computer program results from MAIN 1 confirms this hypothesis. 
Adding a damping term to the ship roll-tank coupling terms, 
however, would have a desirable effect at any frequency. An 
additional tank moment caused by the pressure differential 
would be induced opposing the roll velocity of the ship. 
Thus, the controller effects may be included in the linear- 
ized equations of motion for the type of controller described 
by substituting this reduced tank equation into the non- 


dimensional equations of motion. 





2 
à D 
2 DES t n 
6. e cogn EUR Y 
Hyt 
2 ጻ 1) 
+A, [2 + E ቶ 318 = 0 
me H 


The value of 90 is determined by power considerations. 
The energy dissipation, the product of assumed pump pressure 
and tank volume flow, can be calculated and compared to the 
desired pump power in an iterative fashion until the two 
balance. The effect of pump lag time may be included by 


Seaucting the fraction of the period of roll that was spent 
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BEN snging direction. This factor is about 2mT7/o ~ > fOr 
E Mess than 21/w and zero for E greater than 21/w. 

The above active tank equations are only useful for pre- 
dicting ship motions in a seaway which overpowers the passive 
aspects of the anti-rolling tank. Fortunately, this is the 
case when an active tank would most likely be used. The 
additional damping term contained in : and the need to cal- 
culate mean pump power are small modifications to the com- 
puter program for the equations of motion. Otherwise, the 
same design method as that for the case of passive anti- 
rolling may be used. 

The type of hardware which would give the desired pump 
characteristics is a waterjet type pump situated in the tank 
cross duct. The jet velocity would be much higher than the 
mean fluid velocity, giving rise to a steady pressure inde- 
pendent of average tank fluid velocity. The presure rise is 
caused by the cross duct acting as an inefficient diffuser 
and the pressure reverse may be accomplished by switching the 
waterjet outlet from a nozzle in one direction of flow to an 
opposite nozzle. The actual pump motor would then not have 
to reverse at all, and the lag time would be associated with 


the time required to switch nozzles. 
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VI. Design Examples 

In order to validate the computer program, verify the design 
steps postulated, and investigate the possible range of the solu- 
tions, it was decided to design a series of antirolling tanks 
for a sample family of ships. The ships were chosen to be able 
to illustrate as many aspects and implications of the equations 
of motion as possible within time and effort constraints. 

Three ships were finally settled upon to show basically the 
effects of ship size, degree of sway coupling, and tank size. 
They were a 3600-ton small naval escort type, a 14,500-ton medium 
sized naval vessel which resembles a small aircraft carrier 
or amphibious landing ship, and a 20,400-ton commercial stores 
carrier. The three ships were assumed to be existing designs 
to which an antirolling passive tank was to be added. The two 
naval vessels were possessed of a large degree of sway coupling 
introduced by the parameter Ab = BG/B because they were ships 
with high centers of gravity relative to their other dimensions. 
The commerical stores carrier had only about half as much sway 
coupling, being more typical of a cargo vessel which does not 


possess the high topside payload weights of a naval vessel. 


ie Ship Size 


The principal dimensions of the three vessels were as follows: 


Slip No" I Ship No. 2 SIN = 
^ 3600 14,500 20,400 


L 408 585 500 
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Figure 10B Midship Section Ship No. 2, With Tank 


Locations 
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ከ ት ን I Ship No. 2 Ship No. 3 
B 45 -76.67 71.4 
T 14.5 2]. 3 28.6 
Cb 0.4733 0.5304 0.70 


EE Tank Size 

Based upon the empiricalformula of Section V.2, an approxi- 
mate required value for tank capacity, Ymax., was determined. 
The vertical layout of the ships, by decks, is shown in Figure 10. 
Tank dimensions and parameters were calculated using the collec- 
tion of equations shown on figure 12 by an iterative process. 
The T value for all tanks was maximized in the following way: 
An assumed A, was used for the low tank and applied to the value 
of metacentric height, GM, for the low tank. The residual stabi- 
lity of the reduced metacentric height given by GM (1 - እ,) 
was kept as a required value to be achieved by the other tank 
locations. Since the higher tank locations resulted in a higher 
center of gravity for the ship and therefore lower GM, their 
AL values had to be necessarily lower in order to achieve the 


required residual stability GM (1 - እነ). 


Ship No.1 SINO SN Ship No: 3 


Y max. = .36/logio ለ 5.89 5.09 4.89 
_ 2 

Ymax. = .le thw 5/9 Dez 5.49 ር 

64 (1 - 34) 2.69 5.17 5.25 


Of the three deck locations possible for tanks on Ship No. 1, 


it was found that the tank pluced on the main deck had too small 
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a capacity Y max. and ship-tank coupling A, to be E cal. 
EN the lower tank"and the center tank were further investigated. 

On Ship No. 2 there were only two tank locations possible, 
assuming that keeping the center hangar deck area clear was 
imperative. The higher tank was slightly cramped in order to 
fit it between decks, while the lower tank should probably be 
expanded to meet the hangar deck bottom to eliminate the small 
wasted space between the tank top and the hangar deck. 

On Ship No. 3, the presence of vertical cargo access hatches 
and cargo storage limited again the tank location to two possi- 
bilities: a low tank with the cross duct underneath the cargo 
floor, or to a second deck high tank of small dimensions to 
accommodate the cross duct between cargo hatches. 

For the naval vessels, Ships No. 1 and 2, the tank fluid 
was assumed to be fuel oil such as diesel oil or JP-5 (0? = 1.69 
slugs/ft.?) because such vessels carry a generous reserve fuel 
supply for extended range purposes. In effect, this made the 
passive antirolling tanks a reserve fuel oil tank which is pos- 
sible to ballast with seawater in an emergency. For the com- 
mercial vessel, Ship No. 3, the tank weight was too large to 
FI JC Or rt to be carrying expensive fuel oil, so seawater (P= 
1.99 slugs/ft. >) was assumed to be the tank fluid. This re- 
duced the tank volume somewhat, compared with a fuel filled 
Bank, but ne tank weights of between 1 - 3 percent of the 
ship's total displacement was a factor against their being a 


practical commercial installation. 
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Be Operational Criteria 

The choice of speed and sea state was also made to illus- 
trate the range of solutions of the equations of motion. The 
smallest ship was examined in a sea state 5 TAS = 10 ft.) 
at its operational speed of 25 knots, a typical requirement 
for an escort ship of its class. Ship No. 2 was looked at 
in sea state 6 (u1⁄3 = 15 ft.) at a lower speed of 20 knots 
because it can be expected to behave better than its lighter 
counterpart, and has a lower designed speed. Ship No. 3 was 
examined at its economical cruising speed of 16 knots ina 


1/3,. This is because the tank 


storm of sea state 7 (20 ft. H 
should be able to mitigate motion effects on the cargo in ex- 


Reme conditions. 


4. Initial Program Values 
An initial estimate of the worst heading angle, the point 


for optimizing roll angle versus HL and K was made by takinq 


ተ" 
the point where the peak of the excitation spectral frequency 
distribution corresponded to the higher frequency peak of the 
roll response, assumed about 130 percent of Wo Later work 
Showed that this was not consistently a very good first guess, 
and that a better first estimate could be gotten by running 
Main 2 program with Kt = -0.8 and taking the worst heading angle 
directly from the results. 

The initial values for heading angle and other parameters 


for the first iteration were as follows for each of the ship- 


tank combinations. 


BG 
GM 


At 
Ut 


2 
ust 


Kp 
| Y /m] 
[Ke /Ix| 
D1 
D2 
D3 


Fl 


143 


KG 


Ko 


Ao 


Ao/A, 


LO/R 


Ship No. 1 
low center 
° ISI 
33 6 32.19 
222192. 526 7 

221) 1,55 
1201 OL 
4.11 4.10 
2E NM 

. 67 67 

220 26 

-8250 -8250 
-21450 -21450 
— 192 0097 270500 

0 0 
. 684 .684 

1 1) 1205 

90° 90° 

43 43 

10 10 

1 ከሙ ስ ር 
-2.71 -2.57 
E5213. .5271 
]24.5 91.3 

5.0 5.26 

20 20 

. 498 Ng 


Shap No. 2 
low high 
16.5 IE 
7.39 6.67 
8.24 8.24 
.30 205 
202 . 97 
A572) =24 61 
E NODE 00: 
. 53 . 53 
.41 .41 
-30200 -30200 
-78500 -78500 
-37700 -37700 
0 0 
.684 . 684 
1.05 1.05 
118° eS 
34 34 
15 15 
30.4 511 
-24.0 -21.66 
.4456 .4233 
607 357 
4.61 4.85 
33 35 
.303 .286 


Ship NO. 3 
low high 
| 06፡15 
የከ s 
E tas 
.25 0.075 
1,0 9 
M SES 
A 500 (x 105 
.88 .88 
TE LE 
0 0 
-96790 -96790 
-26800 -26800 
0 0 
8 ከጠጠር 
JG qus 
509 509 
27 27 
20 20 
A 2 06 
gg ae 
.505 . ,493 
610 214 
M NES. - 
> 32 
nos 
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SHIP: TANK: | 


iteration litera tion Ha terna. 








(L/R + A,/A,) 


2240 


KG' = (AKG + 6W6z)/A 


GM” = KG + GM — KG 





Figure 12. Tank Sizing Equations - Worksheet 
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Ship No. 1 Ship No. 2 Ship No. 3 

low center cn high low high 
Weight 8 4 64.6 506 299 611 163 (tons) 
Ze /R . 343 -.339 2436 -. 78 223 zen 
Wr 2097 1.067 900 0602. 1 15,04] 14.75 


5. Summary of Results 


^ 


The results of calculations for " as a function of Wy at 


K. (best) are shown on figure 13 for the first three iterative 
steps of looking at an assumed point in the sea state - heading 
angle space. A discussion of results for each of the three 
ships follows this. It should be repeated that the intent was 
not to design a particular antirolling tank or to apply any 
arbitrary criteria, but instead to get as broad a look as pos- 
Sible at the range and implications of the equations of motion, 
and the solution method. 

Because the center tank of Ship No. 1 did not seem to offer 
any better performance than the low tank, it was not carried 
beyond the first iteration. This low tank in Ship No. 1 became 
the most difficult of the tanks to design because the plot of 


¿1/3 


surprises than the other tanks. At each iteration, the sub- 


versus M at K, (best) exhibits more curvature and other 


E 


optimized parameters which worked well at that particular head- 


ing angle were poor at different headings. Finally, was 


Hr 


chosen as 1.01 by inspection of all the data, and Ky was deter- 


mined by running Main 2 program for short crested seas through 


K. =#0.4, —-0.8, and -1.0. The combination E E O Ky = -0.8, 
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which was the first iteration guess, was just about optimum. 
B crablt performance of the antirolling tank compared to the 
ship without a tank was not spectacular, but because the tank 
uses reserve fuel oil as a working fluid and is partially sit- 
uated in normal tankage space, it was a reasonably successful 
design. 

On Ship No. 2 the low and high tanks were both carried 
through to a final iterative step. The initial guess of ie 
as a worst heading angle delayed the solution somewhat because 
at that heading the characteristic of 9173 did not flatten out 
as expected and also 118° was nowhere close to the real worst 
heading. Both high and low tanks performed well and it is 
interesting to note that in spite of the much smaller size, 
weight, and AL of the high tank it performed as well as the low 
tank. This was probably because the parameter m ror the 
Moni tank had a negative Sign, indicating that roll acceleration 
puached the tank fluid in the direction to give a positive tank 


righting moment. For Us OIN sid eNsicn, the case for 


de 
any tank below the ship center of gravity, the opposite reac- 
ao: tank fluid to roll acceleration occurs. 

On Ship No. 3 only the low tank was investigated fully, 
but it can be noticed at the 50° heading angle that the high 


tank's performance is nearly as good as that of the larger low 


2 
Fr 


A relative comparison of the predicted solutions of the most 


tank. Again, this was due to the negative sign on He 


important tank design parameters is in order. 
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chip No. l Ship No. 2 SO E Ship ፐን 


low tank low tank. high tank low tank 
Y max 5.7° 5.29 ~ E 
እ, 0.20 0.30 0.19 0:25 
M 7201 1 02 25520 12201 
Hour 2303 2.72 15. 50 2.06 
Weight/100A 235 3.49 1457 3.00 
K. =0 8 -1.0 -1.0 -.4 
€ DS 0.26 0.25 0.26 


1/3 UID max (no tank)) 


A performance parameter g (£ =Q 
ELM ed at the bottom of the list, but it is difficult to attach 
much relative meaning to the parameter because all of the ships 
had different operational criteria. Even if the criteria were 
the same, the size differential in the ship's displacements would 
have made comparisons meaningless except in the case of the tanks 
on the same ship. The tank capacity, Ymax, did not enter into 
the optimization because all of the significant tank angles cal- 
culated were less than their saturation angles, the arctan (L/R). 


Of the parameters given above, the first four, Ymax, Has À and 


1... 
Hst are the most fundamentally important because they determine 


the antirolling tank size, dimensions, weight, placement in the 
ship vertically, and they cannot be changed once the installation 
is completed. Ky however may be readily adjusted to suit the 


Particular ship and the other two parameters relate to performance 


or economic matters. 





Be Analysis of Results 
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a. Tank Capacity. The question of how well the tank capa- 


City was utilized may be answered by inspection of the results 


of significant tank angle at the worst heading in short-crested 


seas. The predicted tank angle, 


nt>, may be compared to the 


tank saturation angle, arctan (Lo/R), to give a measure of how 


much capacity was used at the design conditions. 


Statistics were looked at. 


Two tank angle 


is a good engineering estimate 


of the point at which saturation effects must be taken into 


account in calculations. 


was arbitrarily taken as an 


estimate of when tank saturation effects might just start to be 


noticed by the crew of the vessel because of the slapping of the 


sank fluid against the top of the tank. 
17 1919 


the MAIN 2 program and n 


random process to be 3.34/2.00 times nm 


Sma p No. 
low tank 
e | 9.819 
' 16.38° 
arctan (LO/R) . 26.47” 
2 0.47 
p 109, rctan(Lo/R) 0.62 


Ship No. 
low tank 
poy 
11,812 


"OEC 


1/3 


Sho Noc 
high tank 
242 
12.099 


SG 


Ts an Output or 


is found by the statistics of the 


Ship No: 


low tank 
13.489 
ጋጋ ከ 
1997270 
0.89 


1፡1. 


An immediate observation is that the low tank on Ship No. 3 


is the only one which is making good use of its tank capacity at 


design conditions. The other designs are all wasting space by 


building too large a tank for the requirements stated. 


3 
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DE ላ. The ship-tank coupling parameter did not by itself 
seem to ftis results much except as it contributed to the 
tank capacity. Ship size and tank placement had a more definite 
effect. 


er The tank tuning parameter had a secondary effect 


Hy 
on the results because a higher He corresponding to a higher Oy 
is realized by decreasing the area of the cross duct, Al. The 
lighter cross duct weight allowed for more of the fluid to be 
in the wing tanks. Consequently, a change in tank tuning to a 


higher tank frequency made for a tank of slightly greater capacity. 


els Hope The ship~tank decoupling frequency parameter was 


eM t for the reasons previously discussed concerning it. 
Amore negative value of JUNO which is accomplished by moving 
the tank higher vertically, is always associated with better tank 
performance for a tank of a given À. and i capacity Ymax. 

e. Weight Fraction. This parameter can be directly related 
to the cost of the tank. A complementary parameter, the space 
fraction, was not given because the data was unavailable for 
useful ship volume. 


f. K... The tank damping parameter falls within the narrow 


range of sieut -0.4 to -1.4 for all successful designs. This 
Statement is based on analysis of much more data than is presented 
ene” results section. A convincing reason for the truth of this 
s C Int can be found by considering the critical damping for 

a Single degree of freedom system. For a single degree of freedom 


system critical damping implies that there is no overshoot and 


no oscillatory response to a step input. Such critical damping 
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has the Ky value of -2.0. Therefore the ratio of tank damping 
to Bee damping falling in thefrange of 0.2 ter 7 is rea- 
sonable if the tank is to act as a dynamic vibration absorber. A 
K. prediction of greater than -2.0 tells us that the given tank 
parameters are not consistent with a stabilizing tank. An 
example of Ky greater thank ?.0 aseshowamomebigure Vomror the 
fowetank on Ship no. l. 

g. Tank Placement. There is another aspect of vertical 
tank placement in addition to comments already given about how 
it affects hye Ue and GM. A special case of tank operation is 
whenever the duct is closed for some reason while the tank is 
full, for example for emergency maintenance. On high tanks 
especially there is a large GM loss due to the KG rise associated 
with the large tank fluid weight high in the ship. Figure 16A 
gives for an example the high tank on Ship No. 2 to show the 
undesirable effects of this operational condition. 

h. Ship Size. It seems obvious from the results to state 
that a lighter displacement ship such as Ship No. 1] apparently 
is not as well suited a candidate for an antirolling tank instal- 
lation as a larger ship. Whether this is due to the center of 


buoyancy being less submerged or to the inertia of the ship 


allowing for a higher natural frequency remains to be answered. 


7. General Conclusions 

Although none of the design examples is complete in itself - 
primarily because a range of speeds and sea states was not looked 
at for each ship - it was possible to draw a number of general 


conclusions from them. 
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First of all, the equations of motion seem to be valid. 
Uncoupling sway by setting Ab = 0 and res = 0 gives the clas- 
Sical fourth order differential equation for a coupled ship-tank 
system. The assumed excitations are the primary departure from 
the latter, along with the sway coordinate, and hence introduce 
the greatest possible source of error because the reactive force 
and moment terms were also taken from classical theory. A lot 
of numerical accuracy hinges upon the two critical assumptions 
Semuncoupling of roll and sway and application of the excita- 
tions at the center of (static) buoyancy. Model testing of some 
sort is needed to prove or disprove these conjectures. 

The general design method seems to be flexible and effi- 
cient within the range of interest. This was due to the veri- 


Aisa tion that 91/3 versus u, and K a 'drawable' three dimen- 


27 
sional function, is a relatively flat function. This was due 
Es being an integral function. 

The solution of the fifth order differential equations of 
motion by the design method predicts tank parameter values which 
are quite similar and in some cases identical to the classical 
fourth order approach. The roll responses of the two approaches 
is not significantly different. However, a more refined solu- 
tion is possible by the methods of this paper. For example, the 
huuhsrank of Ship No. 2 with ML = 1.2 had a lower weight than 


that for the classically designed tank of Hy = 0.97 due to the 


smaller cross duct associated with a higher tank frequency. 





88 


VII. Conclusions and Recommendations 


1. Conclusions 

The design method proposed in this paper seems to offer an 
attractive, efficient and flexible method for the design of pas- 
Sive antirolling tanks. The analysis for active antirolling tanks 
led instead to a method for a feasibility study for the use of 
active tanks. The active tank design method predicts perform- 
ance of an economical active tank and will select optimum tank 
parameters for it. However, the design of the pump and the pump 
controller was left as a separate step. The active tank analysis 
also demonstrated why active tanks are such unpopular creatures. 
The small difference in tank angle which the pump is able to 
achieve and the small amount of active damping which it adds to 
the equations of motion are not nearly enough to overcome the 
chief cause of tank efficiency suffering in high sea states -- 
the saturation of the tank which places an upper limit on the 
tank righting moment. 

The design method for passive antirolling tanks is attrac- 
tive because it uses a good representation of the real ocean 
wave process in predicting roll angle. Since there is much 
available information on ocean conditions for specific ship- 
ping routes or operational areas, this information may be sub- 
stituted for the assumed Pierson-Moskowitz function. The design 
method is efficient because the equations of motion lend them- 
selves to a complex algebraic form of solution, thereby resulting 


in a high numerical efficiency for the computer program. Because 
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0 01 a large range of solutions may be investigated at little 
cost. Finally, the design method is flexible inithat tankide- 
sign may be approached from a wide variety of given requirements. 
For example, meeting some criteria for operational use, meeting 
Space or weight constraints, or any combination of these may be 
imposed as design requirements. Also, the equations for the 
hydrodynamic and tank coefficients may be changed as necessary 
in the main computer program without altering the subroutine 
which calculates roll response. 

The prediction of roll angles seems reasonably well in agree- 
ment with all of the references, both in magnitude and in dis- 
tribution over heading angles. The short-crested sea predictions 
should be more accurate than those for long-crested or regular 
waves because linearity is better preserved. Only small compo- 
nents of roll for each heading angle are added together by super- 
position, including the worst heading at resonance peaks. 

A general conclusion concerning design parameters is that 
a reasonable passive tank design for preliminary design purposes 
from Appendix B, 


may be estimated for a given À or As = 0.20, 


e Hst 


from the Chadwick-Klotter equation on page 53, K, = -0.80, 


Hr t 


and the tank dimensions from Appendix B. 
2. Recommendations 

It was decided to first investigate and validate the equa- 
tions of motion by taking the simplest case of a passive tank 
and linearized coefficients. As it became apparent that the re- 


sults of this approach were promising, a design method for passive 
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tanks was devised. A thorough checking of the results by model 
or full scale testing seems to be in order as the next step. 

The numerical integration technique seems to offer the best 
meee FOr modelling a first order, or higher order control effector 
EmNCongunction with the fifth order ship and tank equations. 
However, as previously mentioned this technique is inefficient 
in that there are so many integration steps to be performed in 
order to do an optimization, and several orders of magnitude 
more integration steps to be able to look at short crested seas 
Over a range of sea states and ship speeds. In short, the 
numerical analysis technique has a certain merit for analysis 
but 1S a slumsy synthesis tool. Now that the basic equations 
rime reactive forces and moments have been validated, there 
may be a way to overcome this difficulty. This is to combine 
the spectral density function and the spreading function before- 
hand to make again an irregular, short-crested sea by taking 
all the harmonic frequencies and superposing them. A five or 
ten minute real time segment of excitations could be made and 
stored as library data. This data could be manipulated for 
Speed, heading angle and sea state to form a simulation of the 
excitations. This could be used to drive the equations of motion 
with the effector included. Data for events such as significant 
F Te could be obtained from the statistics of the results. 
It would be interesting to see how well the results of such a 
simulation compare with the complex algebraic solution. There 
is room for question as to whether the simulation would be as 


accurate as the complex algebraic solution because of the reduced 
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form some of the linear coefficients would assume. For example, 
the sway damping term must lose its frequency dependence in 

this type of approach because the excitation would be a series 
of absolute quantities derived from a frequency spectrum and 
would lack the frequency information necessary for calculation 
purposes. 

It was noticed when comparing the results of a passive tank 
design with the design solution contained in reference 13 that 
there was very close agreement in overall results in short- 
crested seas except for the case of the region of heading angles 
near 0° and near 180°; that is, head seas and stern seas. The 
author believes that the discrepancy most likely could be due 
to the simplistic approach taken towards obtaining a heading 
angle correction in SectionIII.7. This author assumed a constant 
distribution of wave excitations over an effective shiplength. 

It 1s recommended that other distribution functions more closely 
resembling that which is actually found aboard the ship be tested 
to see whether a closer agreement could be realized. The impli- 
ans of using a second or third order parabolic distribution 
would be worth examining. Before doing this, however, it would 
be wise to carry through on the first recommendation given in 
this section and check on which of the predictions of roll angle 
are closer to model or full scale experimental tests, those of 
this author or those of reference 13. 

A final recommendation is that some of the steps for passive 
tank design could be better automated. For example, the inte- 


grals S' and S" could be accurately handled for any situation 
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by the proper numerical program. Also, the burden of solving 
the tank sizing equations on Figure 14 which need an iterative 
convergent solution should be removed from the tank designer. 
This would be no simple task of programming, however, because 
of the great number of combinations of constraints such as Qu 


Or Wo, being limited by tank dimensional constraints. 
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Appendix A Anti-Rolling Tank Geometry Factors 


He was decided to construct a parametric family of Anti- 
Rolling tanks using the simple tank geometry of Figure A-l. In 
practice, a much more complicated geometry could conceiveably 
exist due to the need to use the ship's bulkheads, sides, and 
decks as part of the tank structure. Also, depending upon the 
tank damping required, it may be necessary to fair in sharp 
corners or other flow obstructions. It is expected, however, 
that the assumed geometry will be close enough to any required 
geometry in any specific installation so that equivalent dimen- 
sions and equivalent integrals (S' and S") can be quickly cal- 
culated and overall dimensions easily estimated. 


The basic geometry is shown below. 
| 


dE — s= xw 





Figure A-l.  Anti-Rolling Tank 
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mie tank form factors are S' and S", 
S' = o] &o/a)a1 


Jara 


I! 


s" 
Performing the required integrations by parts over each 
of the regions of constant sectional area, Aor Ai and እዞ and 


dropping the nearly negligible higher order terms gives the 


"ከ ምርር results. 


~ 


2 ጋ. ስ ”ፓ” ስ ፡ .ያ ጉሜ ZEHN TIU ZUM) 


S" 


s" 


[ 


2R(2L /R + z,/R) 


The two tank natural frequencies are 0, and o ,. 


€ st 
Wy = Y29/S' 
Der =V 29/9" 


S 


ሠ. Vg : IO) ERAS Bao 


Ost “ሀ9/ኳ V are, + 2t/n) 


The above steps have reduced the problem of selecting a 

tank size to selecting a mix of the tank parameters R, (A/A) y 
(L ZR), and (z,/R). It may bewseem that themother tank dimen= 
sions are known once these are established. The range of 
parameters for which this analysis is valid is approximately: 

0.2 < Lo ZR < 1.0 

= z,/R s 

R < B/2 


Notice that for z «-2L. the value of o will be imaginary. 


ile st 


Later work with the equations show that this indicates a 90° 


phase shift of the motion induced in the tank by the ship roll 
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motion, but does not imply any dramatic change in the equations 
eo: motion. 

Figure A-2 summarizes the relationship between the tank 
frequencies and the tank parameters. For most combinations of 
these parameters, and in the range of validity previously men- 
is r, Figure A-2 is sufficiently accurate for preliminary 
tank design purposes. 

A convenient way to determine tank size is by fluid weight. 

Tank weight (fluid) z 20A R(L /R 4 A,/A,) 

An effectiveness paraneterof tank stabilizing potential is 
the ratio of tank righting moment to wave upsetting moment per 
unit of tank or waveslope angle. This ratio, hee 1s also equal 
va the percentage reduction in static GM that the tank induces 
on the ship. 

At = 2pgAR°/K¢ = 12A. R^/35AGM (^ in long tons, Y- p/pH^0) 

Another measure of tank stabilizing potential is the maxi- 
mum waveslope that the tank could statically stabilize. Since 
the maximum tank angle possible is about L /R (radians), the 


static tank capacity Y max is given by: 


Y 
Y 


"€ AtL,/R 


21A L R/35AGM 
max 0 0 
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Appendix B. Coefficient Estimates 


I——Dvdrodynamrc Coefficients 

Figures B-1 and B-2 display the hydrodynamic coefficients 
for ellipsoidal bodies of revolution, as theoretically pre- 
dicted by potential flow theory. The strategy in using these 
figures is to assume a body with the same dimensions (L,B,T,) 
as the ship. In cases where the mass of the body of revolution 
is quite a bit different than that of the actual ship, correction 
can be made by taking the average value of the coefficient for 
a body of the same dimensions (L,B,T) and for a body of the 
same mass as the ship with similar geometry (L/B, B/T, L/T). 

ን ጋ ጋ; ኘሣ 

These are both constants. For order of magnitude 

estimates, Yý = -0.9m, and Ke = -0.25 I. 

Ba Ky and LM 


K given as B44 on the figure, is assumed to be a 


p’ 
constant mean value for reasons previously discussed. Boo on 


the figure corresponds to Y, (0,U) by the notation of this 


i 
A it 1s assumed to be well enough represented as a se- 
quence of linear segments. Figure B-3 graphically depicts this 
linear segmented approximation. The coordinates (Fi, D.) are 
used as computer program inputs. The function is assumed 
constant after point 3; although this is certainly grossly 


inaccurate, the roll response operator for the normal range of 


ships is so small in this region that any assumption about the 
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behavior of the actual function will have essentially no effect 
on the roll angle solution. 
The sway damping component Y 9 (0) is estimated by the 
methods of reference (3) to be: 
E E Yao nc amp 2) TU 
he lesser absolute value is associated with a bare hull 
only, and the factor of two is added for a ship with a large 
amount of rudder surface area and other stabilizing underwater 
Sail areas such as skegs, keels, etc. A particular ship may 
be somewhere between. For a normal amount of rudder area it 
will probably be close to the larger absolute value. 
C. Ko 
Kọ may be found directly from ሺ ም Kö = -mg GM. 
2. Ship Parameters 


A. tp: the center of buoyancy 


The center of buoyancy is accurately estimated by 
using Morrish's Formula for KB. 
tp = T — KB 
KB 


| 


17 315/72 DZ PS DEN) 
Awp = Cwp * L °> B 
Cwp = 0.237 + 0.833 Cp reference 14. 
B. GM and G. 
ive [Cal center of Gravity and GM, rf not known, 


may be estimated together by the relations: 
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KG = KB + BM - GM 
BM = Cit L(B)°/35A 
Cit = -0.013 + 0.10 Cp reference 14 


GM = 0.08 B (See reference 3 for a list of repre- 
sentative GMs) 


Notice that if one or the other of G or GM are known, the 
unknown quantity 1s very accurately determined. Otherwise, 
both quantities suffer from being subject to 'guesstimates'. 

E 
The transverse moment of inertia of the mass of the 


ship may be closely estimated by: ከ = 0.38 B; Ix = m (5.)” 
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Appendix C. Computer Program 


1. General Information 

A program which provides a numerical solution to the 
equations of motion was written in standard Fortran language. 
It is compatible with any Fortran capable system as long as 
the input device number KI and the output device number KO are 
correctly specified in the program. A plotting subroutine 
PICTR was used in the program; if this library subroutine is 
unavailable, replace the CALL PICTR cards with dummy CONTINUE 
cards. 

There are two basic programs, Main 1 and Main 2, 
which are separable decks and need not be run concurrently. 

A single subroutine RESP calculates the roll angle per unit 
waveslope, the tank angle per unit waveslope, and the phase 
angle between the tank and roll. A library subroutine, PICTR, 
jms called to give a pictorial display, but because the results 
are printed out anyway, it is unessential to the program. 

Main 1 section of the program computes optimum tank 
Gamping and tuning for a given set of conditions. 

Main 2 section of the program computes roll response 
versus heading angle in short s ፡፡ seas for a given ship, 
tank, and sea state. 

Base there are many untried variations to the 
program it is advised that anyone using it should proceed with 


care and check results at each step of the computations. 





Computer Program Symbols 
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The following symbol list cross references computer 


variables with the symbols used elsewhere in this text, where 


MAIN PROGRAM 


indicated. 

IX Ix 

LB Àb 

KFI Kd 

KPD Kp 

KY Ky 

KS Ko 

LY Ay 

KE K, 

H3 fee 

LE Le 

@ Spreading Function 
Coefficients 

XLAB Axis Labels 

TITLE Title Tables 

XMT ut 

XXMT Dummy yt 

EE 1⁄2 

RAO Response Ampl. 
Operator 


ሀን 


G Or SK 


AB 

E 
TANG 
KO 
KI 
BG 
GM 
CB 


AL 


TB 


AYVD 
ALT 


Yv2 


ጋ | 
DKTL 


XMST2 


W 2/&w 

Dummy holder 
Tank Angle nm 
n 1⁄3 

Output device number 
Input device number 
BG 

GM 


Cb 


Ky initial value 


Ky increment 
mee 





DMT 
Wi 
DW 
XKPD 
BETA 


DI = D3 
< F3 
YV 


YVD 


Gl 


LT 
LO 
LY 


LB 


M2 
MT 2 
MST2 


KY 
KS 


KT 


MAIN PROGRAM (continued) 
ut increment G2 
0 initial value WS 
w increment, Aw WE 
|KB/1x | W 
P x1 
s 
| አሰድ 는 CTB 
yy ETA or E 
Y. YV1 
2 5 

SUBROUTINE RESP 

Pa W 
AO WS 
AY Gl 
Ab G2 
u TOCBL 
y Aij 
Wt ne 
ME 31 ገ 
Ky Ci 
Ks RX 
K RT 
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) 2 


We 
WO 
X 


Heading angle cor- 
rection el/2. x 


2100 "9 


n 


Y 
V 


Number of increments 
of Aw 


WO 
WS 
Y 1 
(2 
Phase angle 


Matrix elements, 
real 


Matrix elements, 
imaginary 


Intermediate matrix 
steps 


Q/U 
n/v 
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8. Computer Program Input 
The following input data are required by the program. 


Paap Dimensions 


L Ship Length 

B Ship Beam 

T | Ship Draft 

Cb Block Coefficient (Immersed volume/LBT) 

BG Distance from center of gravity to center of buoy- 


ancy, tn Z direction 


GM Metacentric Height 
th Draft to the center of buoyancy 
Le Effective shiplength (between 60% beam sections) 
Tank Parameters 
nt Tank tuning (wt/ws) 
Aut Tank tuning increment for search routine 
Kt Tank damping 
ለጁ Tank damping increment for search routine 
ust" Ship-tank decoupling frequency; primarily a function 
of vertical placement 
At Ship-tank coupling parameter -- $ loss in static 
GM of ship due to tank free surface effect 
Bewatıon Coeffieients 
E Roll damping 
|Y,/m| Added mass coefficient in sway. 
|K,/1x | Added inertia coefficient in roll. 
Dl - D3 Sway damping ordinates for Y] (0,0) 


Fl - F3 Sway damping corresponding to D. 
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wl Initial value (lower cut-off frequency, about 0.2) 
Aw Increment of w (notice, 30 increments are used un- 
less N is changed in the program; the 30 increments 
plus 01 will determine the upper cut-off freq.) 
gi ? Significant waveheight of seas. 
B Heading angle of ship with respect to the seas 
(180' for head seas) 
U Ship speed in the x-direction 
Program Control Variables 
There are three control variables which must be specified 
on the input deck. They are used to switch various functions 
in the program, and to terminate the computations. 
(+1, no PICTR) 
KSTOP +1 Read Type 1 data cards(-1, PICTR at ut = utl, 
Kt = Ktl for RAO for MAIN 1) 
(-1, PICTR of short crested 
seas for MAIN 2) 
0 Stop execution of the program!! 
woe. +1, 0 Main 1 Program executed 
EI Main 2 Program executed 
ESTOP +1 Read Type 2 data cards 
0 Stop execution of program!! 
-1 Read Type 1 data cards; switch back to Main 1 
Program 


caused by the above statements. 


Examples of sample data decks will resolve any confusion 


The intent of having the 


Switching variables was so as to be able to run Main 1, Main 2, 


Or both together, 


and at the same time use some of the 





TOS 


subscripted computer variables (matrix elements) as dummy 
holders for intermediate calculation steps. Unless the data 
decks examples are followed, there are probably several non- 


workable combinations of the switching variables. 


Data Card Input Deck Format 

The general format for all the variables is floating point 
decimal. the exceptions are JSTOP, KSTOP, and LSTOP, which 
are fixed point decimal (12). 

There are three different types of data cards. To run 
Main 1, types 0, 1, & 3 are used. To run Main 2, types 0, 1, 
2, & 3 are needed. The data values on types 1 and 2 need not 
be identical; the reason for having type 2 data cards was to 
permit a number of variations of data without retyping a whole 


type 1 deck. 





110 


Card # Column 


AL I9 711 - 20 21 ㅋ 30 3 - 40 41 - 50 51 = 60 .61.= 70. / ㄴ = 8680 





0 
(column 1-2 KSTOP) 
(column 3-4 JSTOP) 
IA BG GM Cb L B T Ec yi^? 
1B K | y /m| Me blank Le U Key AK, 
1c ሸመ utl A 
ut wl Aw IK,/T, | 8 blank 
1D 581 D2 D3 Fl F2 F3 blank blank 


TE /4&——X-axis label, colums 1-40 





Y-axis label, columns 41-80—h/Y 


2A U jl ሪግ We K. (column 41-42 LSTOP) 43-80 blank 
TE /4———Title to be printed below short crested sea results, col l-80—®/ 
2C /@—— axis lable colums 1-40—2vy«e$——y-axis label, columns 41-80——15»» 


3 Bee el ere. 


Combinations of data cards for successful computer runs: 
1. Find optimum T parameters (Main 1) 
O ee (1A - 1F), 3; set KSTOP - O1, JSTOP - OO 
2. Find response versus heading angle (Main 2) 
INE iF), (2a - 2C), (2A - 2C)... Beer... (Dr 2202 


Set KSTOP = Ol or -l; JSTOP = Ol; LSTOP = OL or -1. 
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be Sample Output Data 


SIG WAVE HT = 10.0 HEADING ANGLE = 90.0 DEG 


PE SPEED 


Mo O00 EFFECTIVE SHIP LENGTH = 26070 


SHIP NAT FREQUENCY = 0.5142 SHIP-TANK COUPLING = 0.2000 


GM= 3.360 SHIP-TANK DECOUPLING FREQ= 1.04 
LB = 0.2167 LY= 8.9732 Gl= 3.3155 G2= 11.5796 
KS = 0.0453 KY=-2.33446 MAS= 231884. IX= 75130336. 
KT FI1/3 i 
-0.400 3.544 
-1.400 5.783 
-0.900 4.870 
-0.400 3.544 
0.100 4.966 
-0.150 4.401 
-0.400 3.544 
-0.650 4.089 
-0.525 3 9233 
298200 3.544 
20827 5 3.706 
MT= F11/3- 3.71 ETA1/3- 5.67 
TANK LINEAR  DAMPING - -0.27 
KT FIl1/3 2 
-0.400 3.474 
21100 4.202 
-0.900 3.602 
-0.400 3.474 
0.100 5.162 
-0.150 4.300 
-0.400 3.474 
-0.650 3.374 
-0.900 3.602 
-0.775 30012 
-0.650 3.374 
MT= FIl/3= 3.37 ETA1/3= 6.33 
TANK LINEAR DAMPING = =0.65 
KT FI1/3 3. 
28:00 7.451 
-1.400 6.939 
2254100 6.430 
-3.400 6.190 
-4.400 6.067 
-5.400 5.996 
-6.400 5.952 
25100 5.922 
MT- fIl/3- 5.92 ETA1/3= 0 32 
TANK LINEAR DAMPING =  -7.40 
MI(MIN)+ 1.35 FI1/3= -1.8098 
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LONG CRESTED SEAS; SPEED= 34.00 SIG WAVE HT= 15.0) 


MT= KT= LT= 
ሽሽ (0) 1.00 0.20 

WS= LB= KS= 
0.40 ክር ታሪ 207 06 


HEADING ANGLE SIG ROLL ANGLE SIG TANK ANGLE (DEGREES) 


0.0 0.00 0.00 
15.0 DS 4 22 
30.0 4.61 9.94 
45.0 6 92 15.25 
60.0 4.44 13:93 
5.0 2.08 6. 29 
70.0 12229 2.99 

1050 0.87 5: 99 
EZ. 0 0.59 2295 
2.0 02539 1295 
150.0 0.24 1.20 
165.0 (LT 0.58 
o0. 0 0.00 0.00 
ORE CRESTED SEAS 
MT= KT= LT= 
1.20 1.00 0.20 
WS= LB= KS= 
0.40 OZ ን. 
HEADING ANGLE SIG ROLL ANGLE SIG TANK ANGLE 

0.0 23:41 UU OD 
15.0 3.42 7.77 
30.0 3.44 8.23 
45.0 3.40 ፡ 
60.0 3.19 OUR 
757, 0 ES 8.08 
90.0 2.14 ou 

ROS. 0 1.48 3.25 
120.0 093 3./2 
1 5520 9 57 መ ነታ 
530.0 0232 2 
i55. 0 0.26 1]... 
180.0 0.23 1.13 
SHIP # 2 HIGH TANK 


z 


zu 
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Appendix D. Other Solution Methods 


ME Analog Solution 

The equations of motion were put into their proper non- 
dimensional form for electrical analog simulation. However, 
it was considered necessary to further investigate the numer- 
ical solution in order to properly normalize the non-linear 
variables, and especially coefficient expressions. Using the 
results of this paper, an analog simulation is a trivial 
problem if the modeling of the coefficients can be 


accomplished. 


e Numerical Integration 

In order to perform this technique, the three couri 
equations forming a fifth-order system must be algebraically 
broken down into five coupled first order differential equa- 
tions. This was done using the set of dimensionless equa- 
tions in preference to the dimensional set. The five state 
variables chosen to represent the fifth-order system were: 
ዕኔ 0, A and n. One of the excitation terms, the th, 
term, was dropped to make the algebra a little more tractable. 
A parametric study was done in order to justify this step. 
The below figure shows that for the test case investigated 
for a normal” sized ship, the iA term has practically no 
bearing on the excitation for ws È 0.2, which was previously 


shown to be the lower boundary of the excitation frequency 
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EA aoronberest. Therefore, the Sway “exci taevon Term was) 
shortened to ~ anderoll exeitaLrom to IAE This kept 


ee excitation term uni-dimensional. 
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Figure Dl. sway Excitations. 


mane Live Coupled Ordinary differential equations, in 0! 


Matrix form, were found to be, 


ጅን e e aA Ee] 
4| 1 0 0 0 0 $ 
d/dt jv" = 12 )/Y) ቾ53/7ቕ2- 23472 0 24772 0 25771 치가 
m 51/53 52/93 52/53 JUL 55/2 n 
ni^ 10 0 0 l 0 n 
8-/፤ 8/8 NEU NS. 
EC "o Eo @ 6 R O 
0 0 0 0 
p|p| 
Z Yo m: 27/8, Bo 26/73 20/72 : 43 
56/* 57/53 59/*3 89/0» 
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The above matrix equation is separated into a ship-tank 
Eunear matrix, and excitation matrix, a controller matrix, 
and a ship-tank non-linear matrix. The coefficients in the 
matrices are all constants for a given set of conditions. A 
Separate equation for p* should be added for a particular 
contoller. 


The identities for the 29 new coefficients, B7B71 21”=0። 





$1799; and 0,707 are given below. 
No እ እ 1 o 
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Figure D2. 
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equations may be drawn as a signal flow graph 
Il os nal Liew grapes useful lor Inter; 
what the basic relationships are between ail 
This information cannot be extracted from tche 


the three coupled equations. 





Ship Anti-Rolling Tank System 
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Simulation of roll by the use of the above equations 
Without the necessity of directly using the spectral density 
function technique could be accomplished by taking an excita- 
tion made of randomly phased harmonic waves, in their proper 
proportions, and integrating for a ship time duration of 
enough cycles, say a thousand, to get a meaningful sample for 
determining es In the interests of amount of effort 
required, this should be done subsequently to finding the 
approximate solution by the complex algebraic method adopted 


by the author. 








thesZ365 
nti-rolling tank design. 





